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Abstract. 

We study the stellar and wind properties of a sample of Galactic O dwarfs to track the conditions under which 
weak winds (i.e mass loss rates lower than ~ 10~** M0 yr~^) appear. The sample is composed of low and high 
luminosity dwarfs including Vz stars and stars known to display qualitatively weak winds. Atmosphere models 
including non-LTE treatment, spherical expansion and line blanketing are computed with the code CMFGEN 
(Hillier & Miller IT998II . Both UV and Ha lines are used to derive wind properties while optical H and He lines 
give the stellar parameters. We find that the stars of our sample are usually 1 to 4 Myr old. Mass loss rates of 
all stars are found to be lower than expected from the hydrodynamical predictions of Vink et al. 1200111 . For 
stars with log ^ 5.2, the reduction is by less than a factor 5 and is mainly due to the inclusion of clumping 

in the models. For stars with log ^ 5.2 the reduction can be as high as a factor 100. The inclusion of X-ray 
emission (possibly due to magnetic mechanisms) in models with low density is crucial to derive accurate mass loss 
rates from UV lines, while it is found to be unimportant for high density winds. The modified wind momentum 
- luminosity relation shows a significant change of slope around this transition luminosity. Terminal velocities of 
low luminosity stars are also found to be low. Both mass loss rates and terminal velocities of low L stars are 
consistent with a reduced line force parameter a. However, the physical reason for such a reduction is still not 
clear although the finding of weak winds in Galactic stars excludes the role of a reduced metallicity. There may 
be a link between an early evolutionary state and a weak wind, but this has to be confirmed by further studies of 
Vz stars. X-rays, through the change in the ionisation structure they imply, may be at the origin of a reduction 
of the radiative acceleration, leading to lower mass loss rates. A better understanding of the origin of X-rays is of 
crucial importance for the study of the physics of weak winds 
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1. Introduction 

Massive stars are known to develop winds so intense that 
mass loss rate turns out to be the main factor governing 
their evolution (e.g. Chiosi & Maeder'1986'). The mecha- 
nism responsible for such strong outflows was first pointed 
out by Milne H1926|l when observations of winds were not 
yet available: the radiative acceleration in these bright ob- 
jects was suspected to be large enough to overtake gravita- 
tional acceleration, creating expanding atmospheres. The 
first quantitative description of this process was given by 
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Lucy & Solomon H1971|l who computed mass loss rates 
due to radiative acceleration through strong UV resonance 
lines. Castor, Abbott & Klein (|1975|l made a significant 
improvement in the understanding of winds of massive 
stars in their detailed calculation of radiative accelera- 
tion including an ensemble of lines by means of their 
now famous formalism and found mass loss rates ~ 100 
times larger than Lucy & Solomon H1971|l . The theory 
of radiation driven winds developed by Castor, Abbott & 
Klein was further improved by Pauldrach et al. : 1 9>( i 1 and 
Kudritzki et al. I|1989|l who included the effect of the finite 
size of the star in the radiative acceleration. 

In parallel to theoretical studies, observational con- 
straints on the wind properties of massive stars were ob- 
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tained. Most methods relied on either the measurement 
of infrared and radio excess emitted in the wind of such 
stars (Howarth & Prinja 1989, Leitherer ll988[ Lamers & 
Leitherer ll993|l . or on the analysis of UV and optical emis- 
sion or P-Cygni lines (e.g. Leitherer 119881 Haser 119951 
Puis et al. 1996). The results confirmed the prediction of 
the theory that the mass loss rate should scale mainly as 
a power law of luminosity (e.g. Howarth & Prinia [T989|l 
and that the terminal velocities are directly proportional 
to escape velocities (e.g. Lamers et al. ll995ll . Another suc- 
cess of the radiation driven wind theory came from the 
so called modified wind momentum - luminosity relation 
(hereafter WLR). Kudritzki ct al. J1995 ) showed that the 
quantity Mvo^^/R (with M the mass loss rate, Voo the 
terminal velocity and R the stellar radius) should depend 
only on luminosity (contrary to M which also depends 
slightly on the star mass) which was soon confirmed by the 
spectroscopic analysis of O and B stars (Puis et al. 119961 
Kudritzki et fd. HM^ . This finding was quite exciting since 
once calibrated, the WLR could be used as a distance indi- 
cator up to several Mpc (Kudritzki 1998). Recent determi- 
nations of wind parameters with sophisticated atmosphere 
codes confirm the good agreement between observational 
constraints and theoretical predictions for bright O stars, 
both in term of mass loss rate (for which the most recent 
predictions arc those of Vink et al. IMHH IMIT|l and WLR 
(see Herrero, Puis & Naiarro OUO^ Crowther et al. 12002 
Repolust et al OUMjl . 

In spite of these encouraging results, the behaviour of 
the wind properties of O stars with relatively low lumi- 
nosity seems to be a little more complicated. Martins et 
al. Ij2002bl 120041 hereafter paper I) have shown that the 
stellar components of the star forming region N81 of the 
SMC are O dwarfs with low luminosities and surprisingly 
weak winds: the mass loss rates are lower than 10 ~* M0 
yr""'^ and the modified wind momenta are nearly 2 orders 
of magnitude lower than expected from the WLR obtained 
for bright stars. Bouret et al. H2(JU3fl also found low mass 
loss rates for the faintest of the NGC 346 dwarfs they 
analysed. Although all stars were in the SMC, we showed 
in paper I that metallicity may not be the only factor re- 
sponsible for such a strong reduction of the wind strength. 
In particular, we showed that a Galactic star - 10 Lac - 
displayed a similar weak wind. One of the explanations 
we highlighted was a possible link with the youth of the 
stars since most of them were (or were suspected to be) Vz 
stars, i.e. young stars lying close to the ZAMS (Walborn & 
Parker|1992). Another possibility was a break down of the 
scaling relations (especially the WLR) at low luminosity. 
This reduction of the wind strength at low luminosities 
was in fact already mentioned by Chlebowski & Garmany 
l|1991|l more than a decade ago. 

In this paper, we try to investigate more deeply the 
wind properties of low luminosity Galactic stars. The aim 
is 1) to see if one can exclude the effect of metallicity to 
explain the weakness of the winds, 2) to test the hypoth- 
esis of the link between the weakness of the wind and the 
youth of the stars and 3) to quantify the wind proper- 



ties of faint O stars and the luminosity below which such 
weak winds are observed. For this, we study a sample of O 
dwarfs with both low and high luminosities. Stars known 
to display qualitatively weak winds are included together 
with stars belonging to the Vz subclass. We selected stars 
showing weak UV lines usually sensitive to winds (from 
the lUE atlas of Walborn et al. I1985|l and/or with low 
mass loss rates from the study of Chlebowski & Garmany 
1)1991(1 . We also included Vz stars (N. Walborn, private 
communication) and bright stars (two in common with 
the Repolust et al. 120041 sample) to examine the depen- 
dence of the wind properties on luminosity. Finally, stars 
from the young star forming region in the Rosette nebula 
were included. 

The remainder of the paper is organised as follows: 
In Sect. 13 we give information about the observational 
data we used; Sect. |31 explains how we derived the stellar 
and wind parameters; Sect. 01 gives the results for individ- 
ual stars; Sect. [3 highlights the importance of X-rays and 
magnetic fields in weak wind stars, while Sect. |0| discusses 
possible sources of uncertainty; the results are discussed 
in Sect. Hand the conclusions are given in Sect.|Hl 

2. Observations 

2.1. Optical 

Various sources have been used for the optical spectra of 
the stars studied here. First, the VLT archive provided 
UVES spectra for HD 152590, HD 38666 and HD 46202. 
The instrumental resolution varies between 0.04 A and 
0.1 A, due to different slit widths. The UVES pipeline was 
used for the reduction of the data. Second, optical data for 
HD 34078 and HD 15629 were retrieved from the La Palma 
archive. Spectra obtained with the instrument ISIS on the 
WHT were reduced using standard procedures under the 
ESO/MIDAS environment. The spectral resolution is 0.9 
A. Third, spectra of HD 93204, HD 93250 (EMMI) and 
HD 15629 (La Palma) were provided by Artemio Herrero 
and Danny Lennon and have a typical resolution of 0.95 
A. Finally for stars HD 93146, HD 93028, HD 46223 and 
HD 42088, we used EMMI spectra obtained during the 
nights of 29, 30 and 31 December 2003 on the ESO/NTT 
in La Silla, under the program 72.D-0038(A) (PI Martins). 
These spectra were obtained in the red mode of the instru- 
ment and provided the Ha profiles. The IRAF package was 
used for the data reduction. For a few stars, we were left 
with several spectra of the same wavelength range. In that 
case, we always chose the spectra with the best resolution. 
The signal to noise ratio depends on the instrument used 
but is usually larger than 100 in most lines of interest. 

2.2. UV 

The lUE archive was used to retrieve the UV spectra of all 
the stars of this study. Spectra in the range 1150-2000 A 
obtained with the Short Wavelength Prime (SWP) camera 
were selected. The typical instrumental resolution is 0.2 A 
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and a S/N ratio of the order of 10. The normalisation was 
made "by eye" and turned out to be somewhat uncertain 
below 1200 A. 

We also retrieved FUSE spectra when available from 
the MAST archive. The data are provided already re- 
duced (without binning) and co-added by the CALFUSE 
pipeline (version 1.8.7), and we simply normalised them 
by eye. Due to the strong Galactic interstellar absorp- 
tion, many broad absorption bands form H2 render the 
bluest part of the FUSE spectra useless for our purpose 
(e.g. Pellerin et al. lMT^ . We mostly used the 1100-1180 
A range which has a better signal to noise ratio than the 
lUE spectra for such wavelengths and extends to shorter 
wavelengths. 

3. Method 

Our main concern is to derive wind parameters (mass loss 
rates, terminal velocities) and modified wind momenta. 
However, such determinations require reliable stellar pa- 
rameters, especially effective temperatures. Indeed, any 
uncertainty in Tcff can lead to an error on M . We thus first 
estimate the stellar parameters using the optical spectra, 
and then we use the UV range + Ha line to determine the 
wind properties. 

3.1. Stellar parameters 

The main stellar parameters have been determined from 
blue optical spectra. As such spectra contain diagnostic 
lines which are formed just above the photosphere and are 
not affected by winds, plane-parallel models can be used 
for a preliminary analysis. Hence, we have taken advan- 
tage of the recent grid of TLUSTY spectra (OSTAR2002, 
Lanz & Hubenv l2()()2|l . This grid covers the log g - T^s 
plane for O stars and includes optical synthetic spectra 
computed with a turbulent velocity of 10 km s~^. The 
models include the main ingredients of the modelling of O 
star atmospheres (especially non-LTE treatment and line- 
blanketing) except that they do not take the wind into 
account (see Hubeny & Lanz I19il^ for details). 
Our method has been the following: 

- V sini : we adopted the rotational velocities from 
the literature (mostly Penny I1996|l and refined them in 
the fitting process when possible. 

- Tcff .- the ratio of He i A4471 to He 11 A4542 equivalent 
widths gave the spectral type which was used to estimate 
Toff from the Teff-scale of Martins et al. I|2002a|l . Then, 
TLUSTY spectra with effective temperatures bracketing 
this value were convolved to take into account the rota- 
tional velocity and instrumental resolution, and the result- 
ing spectra were compared to the observed profiles of the 
He I A4471 and He 11 A4542 lines. The best fit led to the 
constraint on Tcff. As the OSTAR2002 grid has a relatively 
coarse sampling (2500 K steps), we have often interpolated 
line profiles of intermediate temperatures. A simple linear 



interpolation was used. For the stars for which the He i 
A4471 and He 11 A4542 lines were not available, we used 
He I A5876 and He 11 A5412 as the main indicators. 

Secondary Toff diagnostic lines such as He i A4388, He i 
A4713, He i A4920, He i A4144, He i A5016 and He 11 A4200 
were also used to refine the determination (when avail- 
able). The uncertainty in Toff depends on the resolution 
of the spectra and on the rotational broadening. Indeed, 
the broader the profile, the lower the precision of the fit of 
the line. The typical error on Tefi is usually of ± 2000 K 
but can be reduced when many optical He lines are avail- 
able. Note that the errors we give are 2a errors (we have 
Tcff — error < T^s < Tcb + error). 

We also checked that our final models including winds 
computed with CMFGEN fitted correctly the optical 
lines. It turns out that the agreement between TLUSTY 
and CMFGEN is very good as already noticed in previous 
studies (e.g. Bouret et al. "20(13). The problem recently 
highlighted by Puis et al. (,2005J concerning the weakness 
of the He I singlet lines between 35000 and 40000 K is 
in fact related to subtle line blanketing effects and is 
solved when both the turbulent velocity is reduced (in 
the computation of the atmospheric structure) and other 
species (Neon, Argon, Calcium and Nickel) are added 
into the models (see Sect. I4.5| . 

- log g : Fits of the wings of H7 led to constraints on 
logg. Once again, interpolations between the OSTAR2002 
spectra were made to improve the determination as the 
step size of the OSTAR2002 grid is 0.25 in log g. H/3, 
which behaves similarly to H7, was used as a secondary 
indicator. The typical uncertainty in log g is 0.1 dex. 

Once obtained, these parameters have been used to 
derive L, R and M: 

- Luminosity : with Tcff known, we have estimated a 
bolometric correction according to 

BC(Teff) = 27.66 - 6.84 x logT^fi (1) 

which has been established by Vacca et al. I|1996|l . Visual 
magnitudes together with estimates of the reddening and 
the distance modulus of the star have then lead to My 
and L from: 

log^ = -0.4(Mv + SC-Mq) (2) 

the error on Tcff leads to a typical error of 0.2 dex on 
BC. Note that we have recently revised the calibration of 
bolometric corrections as a function of Tcff (see Martins 
et al. ,2005 ) . but it turns out that due to line-blanketing 
effects, BCs are reduced by only 0.1 dex, which translates 
to a reduction of log L by 0.04 dex, which is negligible here 
given the uncertainty in the distance. 

The solar bolometric magnitude was taken as equal 
to 4.75 (Allen |1976|. We want to caution here that for 
most of the stars of this study, the distance is poorly 
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Table 1. Adopted parameters for our program stars. The photomerty comes from Chlebowski & Garmany H1991fl . 
Hiltner H1956|l . Howarth & Prinja H1989|l . Humphreys H1978|l . Mai'z-Apellaniz et al. {200^, Massey et al. H2001|l . 
Nicolet fT^TSll, Puis et al. (1996), Schild et al. ifTUSSjl . Walborn et al. fSOM). Distance modulus are taken from 
DeGioia-Eastwood et al. (^2001,1 . Humphreys l|1978(l . Mai'z-Apellaniz et al. (2004j), Markova et al. H2004fl . Massey et al. 
H2Ul)l|l . Mason et al. H1998|l and the Webda database. References to the observing data sets retrieved from archives 
are also given when available. 
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known (with sometimes a difference of 1 magnitude on 
the distance modulus between existing determinations). 
This leads to an important error on the luminosity. 
As this last parameter is crucial for the calibration of 
the modified wind momentum - relation, we decided to 
take the maximum error on L by adopting the lowest 
(resp. highest) luminosity (derived from the lowest -resp. 
highest- extinction, distance modulus and bolometric 
correction) as the boundaries to the range of possible 
luminosities. The typical error on L is ~ ± 0.25 dex, and 
the main source of uncertainty is the distance. 

- Radius : Once Tcff and L are known, R is simply 
derived from 

where a is the Stefan Boltzmann constant. Standard errors 
have been derived according to 

Alogi? = 0.5V(AlogL)2 + (4Alogreff)2 (4) 

- M : The (spectroscopic) mass is derived from g and 
R according to 



and the standard error is given by 

A log M = V(Alogg)2 + (2Alogi?)2 (6) 

With this set of stellar parameters, we have run mod- 
els including winds to derive the mass loss rate and the 
terminal velocity (see next section). The stellar parame- 
ters giving the best agreement between observations and 
models with winds were adopted as the final stellar pa- 
rameters. 



3.2. Wind parameters 

UV (and FUV when available) spectra and Ha profiles 
were used to constrain the wind parameters. In the case 
where mass loss rates were low, priority was given to UV 
indicators since Ha becomes insensitive to M: for such 
situations, we checked that the Ha line given by our mod- 
els with M estimated from UV was consistent with the 
observed line. We want to stress here that it is only be- 
cause metals are now included in a reliable way in new 
generation atmosphere models that such a study is pos- 
sible. Indeed, UV metallic lines now correctly reproduced 
allow to push the limits of mass loss determination below 
- 10-^ Mq yr-i. 

Models including stellar winds were computed with 
the code CMFGEN (HiUier & Miller HHMIl ■ This code al- 
lows for a non-LTE treatment of the radiative transfer 
and statistical equilibrium equations in spherical geometry 
and includes line blanketing effects through a super-level 
approach. The temperature structure is computed under 
the assumption of radiative equilibrium ^. At present, 
CMFGEN does not include self-consistently the hydro- 
dynamics of the wind so that the velocity and density 
structures must be given as input (but hydrodynamical 
quantities computed from the final atmosphere model are 
given as output). In order to be as consistent as possible 
with the optical analysis, we have used TLUSTY struc- 
tures for the photosphere part and we have connected 
them to a classical /3 law {v = Vao{l — ^r)^) representing 
the wind part. We chose /3 = 1.0 as the default value for 
our calculation since it turns out to be representative of 
O dwarfs (e.g. Massa et al. 12113) • The TLUSTY struc- 
tures have been taken from the OSTAR2002 grid or have 
been linearly interpolated from this grid for Tcff not in- 



Note that in some models adiabatic cooling was also in- 
cluded, see Sect.0 
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eluded in OSTAR2002. This method has also been used 
by Bouret et al. (20(S|) and has shown good consistency 
between CMFGEN and TLUSTY photospheric spectra. 

Clumping can be included in the wind models by 
means of a volume filling factor / parameterised as fol- 
lows: / — foo + (1 — /oo)e "'"it where foo is the value of 
/ at the top of the atmosphere and Uinit is the velocity at 
which clumping appears. As in Bouret et al. H2UU3|I . we 
chose Winit — 30 km s^^. 

A depth independent microturbulent velocity can be 
included in the computation of the atmospheric struc- 
ture (i.e. temperature structure -I- population of individ- 
ual levels). We chose a value of 20 km s~^ as the de- 
fault value in our computations. Several tests (Martins et 
al. I2002al Bouret et al. 12003(1 indicate that a reasonable 
change of this parameter has little effect on the emer- 
gent spectrum, except for some specific lines (see Sect. 
I4.5|l . For the computation of the detailed spectrum re- 
sulting from a formal solution of the radiative transfer 
equation (i.e. with the populations kept fixed), a depth 
dependent microturbulent velocity can be adopted. In 
that case, the microturbulent velocity follows the relation 

^iiirb(^) — ^min ~t- {Vmax ^min) wherC Z^rnm and Vjnax 

are the minimum and maximum microturbulent velocities. 
By default, we chose Wmm = 5 km s~^ in the photosphere, 
and Vmax = min (0.1 Voo, 200) km s"^ at the top of the 
atmosphere. For some stars, we had to increase Vmm from 
5 to 10 km s~^ to be able to fit correctly the observed 
spectra. 

CMFGEN allow the possibility to include X-ray emis- 
sion in the models. In some cases (see next Section), we 
had to include such high energy photons. Practically, as 
X-rays are thought to be emitted by shocks distributed in 
the wind, two parameters are adopted to take them into 
account: one is a shock temperature (chosen to be 3 x 10^ 
K since it is typical of high energy photons in O type stars, 
e.g. Schulz et al. 1200^ Cohen et al. l2003|l to set the wave- 
length of maximum emission, and the other is a volume 
filling factor which is used to set the level of emission. With 
this formalism. X-ray sources are distributed throughout 
the atmosphere and the emissivities are taken from ta- 
bles computed by a Raymond & Smith code (Raymond & 
Smith HUTTjl . We include X-rays in the models for the four 
faintest stars and in test models for the strong lined star 
HD 93250 as explained in Sect. using measured X-ray 
fluxes or a canonical value of Lx/^boi = —7.0. 

The main wind parameters we have determined are 
the mass loss rate (M) and the terminal velocity (vqo). 
Constraints on the amount of clumping were also de- 
rived when possible. The terminal velocities have been 
estimated from the blueward extension of the absorption 
part of UV P-Cygni profiles which occurs up to + 
Vmax where Vmax is the maximum microturbulent veloc- 
ity described above: fits of the UV P-Cygni profiles using 
the above relation for microturbulent velocity allows a di- 
rect determination of Voo- Note that other definitions of 
the terminal velocity exist (see Prinja, Barlow & Howarth 



H1990|l V. The typical uncertainty in our determination of 
Voo is 200 km (depending on the maximum microtur- 
bulent velocity we adopt). 

Fits of strong UV lines such as N v A 1240 C iv 
AA1548,1551 Si iv AA1394,1403 O v A1371 and N iv A1718 
provide constraints on M. Ha is also sensitive to M: in 
dwarfs with weak winds, a quasi photospheric profile is 
expected but the line can be used to estimate upper lim- 
its on the mass loss rate as it is filled by emission when 
M > 10~^ M© yr~^. Also, in the case of weak winds 
C IV AA1548,1551 is actually the only line showing some 
sensitivity to wind and was in several cases our best M 
estimator. Given this, we tried to adjust the mass loss rate 
(and clumping parameters) to get the best fit of both the 
UV wind sensitive lines and Ha. 

As regards the abundances, we have taken as default 
values the solar determinations of Grevesse & Sauval 
fl998) since the stars of this study are are all Galactic 
stars. CNO solar abundances have been recently revised 
downward by Asplund H2()()4|l . However, we preferred to 
rely on the Grevesse & Sauval abundances since they have 
been widely used in previous studies of massive stars and 
are therefore more suited for comparisons. When neces- 
sary, we indicate if these abundances have been changed 
to get a better fit. 

4. Results 

In this Section, we present the results of the analysis for 
each star (from the latest to the earliest type ones) and 
highlight the main difficulty encountered in the fitting pro- 
cess. The observed properties and adopted parameters are 
given in Tabled The derived stellar and wind parameters 
are gathered in Table El while results from previous stud- 
ies of wind properties are given in Table |21 

Spectra from atmosphere models are convolved to in- 
clude the instrumental resolution of the observational data 
and the projected rotational velocity of the star. The 
wavelength range between '-^ 1200 and ~ 1225 A is not 
used in the spectral analysis since it suffers from a strong 
interstellar Lyman absorption. 

A general comment concerning effective temperatures 
is that we find values lower than previous determinations 
(see table 13)) since line-blanketing is included in our mod- 
els. This effect is now well known and has been highlighted 
by several studies (Martins et al. |2 Q02a. 2,005 , Crowther 
et al. 120021 Markova et al. 12004(1 . Our mass loss rates are 
also generally lower than previous determination for rea- 
sons discussed in Sect. 17.2.21 

4.1. HD38666 

HD 38666 (also /i Col) is an 09.5V runaway star for which 
we derive an effective temperature of 33000 K from the fit 
of the optical He lines, as shown in Fig. A value of 
logg 4.0 is derived from the Balmer lines. 

The Ha and UV fits are given in Fig.Eland Fig.El The 
best fits are obtained for M = 10~^'^ M© yr~^ and Voo 
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Table 2. Derived stellar and wind parameters of the Galactic stars studied here. The escape velocities were computed 
from the spectroscopic derived mass, radius and luminosity. The evolutionary masses have been estimated from the 
isochrones of Lejeune & Schaerer (2001). Projected rotational velocities were adopted from Penny fl996T, Howarth & 
Prinja 1)1989(1 . Howarth et al. I|1997|rl and Villamariz et al. l)2UU2f) and refined in the fitting process when possible. 



HD 


38666 


34078 


46202 


93028 


152590 


93146 


42088 


Teff(kK) 

BC 
log g 


33±1 

q nc+O.lO 

-O.ZU_0_Qg 

4.0±0.1 


33±2 

q ncr+O.lO 

-O.ZO_Q Qg 

4.05±0.1 


33±2 

q nr + O.lO 

-O.^U_Q Qg 

4.0±0.1 


34±2 

-3.34l°:i? 

4.0±0.1 


36±1 

q co+O.OS 

-O.O^_0 Qg 

4.10±0.1 


37±2 

-3.59tVil 
4.0±0.1 


38±2 

-3.67l»;l| 

4.0±0.1 


V sim(km s~^) 

R(R0) 

Mspectro (Mq) 
M^^ol (Mq) 


111 

4.66i»;« 
6.58t«;S 

1 (7+0.25 
J^O-0.10 

19 


40 

4.77i»;g 

23t.ll 
20 


30 

8.38tn? 

26t? 

21 


50 

O.UO_Q 22 

q 71+3.11 

' ^-2.38 

q^ + 28 

■J^-16 
25 


66 

4 70+0.33 

^- ' ^-0.24 
(. .2+3.00 
O-*^-2.05 
10+23 

22 


80 

r 99+0.23 
J.^^_0.25 
07+3.29 
».f 1 -2.49 

36tll 
30 


60 

r 9q + 0.19 
J.^<J_0.19 

qq+24 

31 


Vesc (km S~^) 


920 


1043 


1046 


1112 


1015 


1106 


1100 


Uoo(km s"-') 


1200 


800 


1200 


1300 


1750 


2800 


1900 


log M(Mq yr-1) 


-9.5±0.7 


-9.5±0.7 


-8.9±0.7 


-9.0±0.7 


-7.78±0.7 


-7.25±0.7 


-8.0±0.7 


vturb" (km s~-^) 


5-120 


5-80 


5-120 


5-130 


10-175 


5-200 


5-190 


foe 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


log Mvink (M© yr-i) 


-7.41 


-7.38 


-7.23 


-6.97 


-7.15 


-6.58 


-6.17 


log MvooV^ 


24.79 


24.64 


25.44 


25.41 


26.67 


27.50 


26.57 



the first value is the minimum and the second one the maximum microturbulent velocity (See Sect. 13.211 



Table 2. Continued. 



HD 



93204 



15629 



46223 



93250 



Teff(kK) 

BC 

log g 

V sim(km s~^) 
log 7^ 
R (Re) 

Mspectro (Mq) 
Me„o! (Mq) 

Vesc (km s~^) 



40±2 

q 09 + 0.14 

■J-°^-0.16 

4.0±0.1 
130 

c CI +0.25 
^•'-'J^-0.20 

41 

1178 



41±2 

-3.89tl:ll 

3.75±0.1 
90 

O.OD_Q ig 

3otT, ■ 

44 
799 



41.5±2 
qq+0.15 

4.0±0.1 
130 

5.57t°;°^ 



11. 



+22 

16 



10 
+ 1.78 
1.58 



51 

45 

1157 



44±2 

^■J^"-0.12 

4.0±0.1 
110 



+0.25 
-0.17 



6.12 
19.87l^t? 

105 
1461 



Woo(kms ^) 

log M(Mq yr-i) 

Vtur6(km S"""") 
foo 

\ogMvM (Mq yr^^) 

log MVooV^ 



2900 

-6.75±0.7 

5-200 

0.1 

-6.11 
28.05 



2800 

-6.5±0.7 

10-200 

0.1 

-5.74 
28.29 



2800 

-6.5±0.7 

10-200 

0.1 

-5.97 
28.28 



3000 

-6.25±0.7 

10-200 

0.01 

-5.25 

28.68 



= 1200 km s~^. Importantly, X-rays have been included 
in the modelhng with log — —6.87 as indicated by 
the observed X-ray emission (see Table EJ. If this high en- 
ergy component is not included, we need a mass loss rate 
10 times lower to fit the C iv AA1548,1551 fine. The rea- 
son for this is that the ionisation structure of the wind 
is increased when X-rays are present, leading to a lower 
C IV ionisation fraction, and thus requiring a higher mass 
loss rate to reproduce the observed line profile (see Sect. 



O for a more complete discussion) . Note that the fit of the 
C IV AA1548,1551 profile is not perfect. This is due to the 
presence of interstellar absorption which adds to the pho- 
tospheric component. However, the fit of the blueshifted 
wind part of the line is good and is not affected by inter- 
stellar absorption (see also paper I) . Previous estimates of 
M range between lO"®-^^ and lO"'^-^^ Mg yr'^ (see Table 
OJ. Our estimate is lower than all these determinations. 
The determination of Leitherer H1988|l relies only on the 
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Table 3. Comparison between our derived wind param- 
eters (M, Voo, shown in bold text) and previous determi- 
nations. (1) Leitherer (TMS), (2) Bernabeu et al. (1989), 
(3) Howarth & Prinja (1989), (4) Prinja et al. 1^1990,1 . (5) 
Chlebowski & Garmany p991|l . (6) Lamers & Leitherer 
(UnSSIl, (7) Lamers et al. (8) Puis et al. (1996), 

(9) Howarth et al. l)TM7jl . (10) Lamers et al. (|T999,1 . (11) 
Repolust et al. PMji . (12) Markova et al. ^M^ . 



log M 



< -7.22 (1), -7.8 (3) 
-8.31 (5) 

9.5 

-6.6 (3) 
9.5 

< -6.87 (1), -7.2 (3) 
-8.10 (5) 

8.9 

-7.0 (3) 
9.0 

-6.9 (3), -7.36 (5) 
-7.78 
-6.9 (3) 
7.25 

-6.35 (1), -7.0 (3) 
-6.82 (5), -6.42 (12) 

8.0 

-6.1 (3) 

6.75 
-5.89 (11) 

6.5 



-5.75 (1), 
-5.62 (5), 
-5.68 (10) 

6.5 



-5.8 (3) 
-5.85 (6) 



-4.9 (8), -4.6 (3) 
-5.46 (11) 

6.25 



2000 (1), 1000 (3) 

1200 

750 (3) 
800 



2100 
1150 
1200 

1500 
1300 

2150 
1785 
1750 

2975 
2565 
2800 

2550 
2300 
2155 
2200 
1900 

3250 
2890 
2900 

3200 
2810 
2800 

3100 
3100 
2800 
2910 
2800 

3250 
3250 
3230 
3000 



1), 750 (3) 
5), 1590 (2) 



3), 1780 (2) 



3), 2300 (5) 
9) 



3) , 3200 (2) 

4) , 2640 (9) 



1), 2030 (3) 
5), 2420 (2) 
4), 2215 (9) 
12) 



3) , 3180 (2) 

4) , 2900 (7) 



11), 3220 (2) 
9) 



1), 3100 (3) 
5), 2800 (6) 
10), 3140 (2) 
4), 2900 (7) 



3350 (3) 
3470 (2) 



Ha wind emission, which in the case of low mass rates 
is very small and difficult to disentangle from the photo- 
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Fig. 1. Best fit (red dashed line) of the optical spectrum 
(black solid line) of HD 38666. The effective temperature 
is 33000 K, logs = 4.0 and V sini = 110 km s'^. 




Fig. 2 

(black 
10 



9.5 



m 

. Best fit (red dashed line) of the observed Ha line 
solid line) of HD 38666. We have derived M = 
Mq yr^^ and v^o was 1200 km s~^. 



spheric absorption. The studies of Chlebowski & Garmany 

mm 



fl991 ) and Howarth & Prinja (|1989ll are based on the fit of 
UV resonance lines with the following method: the optical 
depth as a function of the velocity (only for unsaturated 
profiles) is determined by profile fitting; from this, the de- 
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Fig. 3. Best fit (red dashed line) of the UV spectrum 
(black solid hnc) of HD 38666. For this model, M = IQ-^-^ 
Mq yr-i, Voc = 1200 km s"^ and logLx/L^oi = -6.87 



Fig. 4. Best fit (red dashed line) of the optical spectrum 
(black sohd line) of HD 34078. Here, T^s = 33000 K, log 5 
= 4.05 and V sini = 40 km s~^ 



termination of the mass loss rate requires the adoption 
of an ionisation structure which may or may not be rep- 
resentative of the real ionisation in the atmosphere. This 
assumption may affect the M determination. 

As for Woo, a higher terminal velocity leads to a too- 
much-extended blueward absorption in C iv AA1548,1551. 
The value of Woo we derive is just above the escape veloc- 
ity. Leitherer H1988|l estimated Voo = 2000 km s^^ while 
Howarth & Prinja ((TMUIl found 1000 km s'^ (see Table 
EJ, illustrating the uncertainty in the exact value of the 
terminal velocity of HD 38666. 

4.2. HD34078 

HD 34078 (also AE Aur) is a runaway 09.5V star pos- 
sibly formed as a binary (with ^ Col, see Hoogerwerf et 
al. I2001|l and ejected after a binary - binary interaction 
with L Ori (see Sect. I4.l|l . Fig. 0] shows the best fit of the 
optical spectrum. From this best fit model, we derive an 
effective temperature of 33000 K. This is confirmed by 
the good fit of the iron lines shown in Fig. Note that 
the presence of C 11 AA6578,6582 in the model (see Fig. 
O may indicate a slightly too low effective temperature. 
Test models reveals that increasing T^s to 34000 K weak- 
ens this doublet. However, since these lines also depends 
on the C abundance, we prefer to rely on the He lines and 
UV iron forests estimate. Also, improving the model atom 
for C II produces a weaker line since over recombination 
routes are available, reducing the populations of the C 11 
AA6578,6582 transition levels. This shows that the value 
for Tcff reported in Table El should be considered with its 




m 

Fig. 5. Best fit (red dashed line) of the observed Ha line 
(black solid fine) of HD 34078. Here, M = lO"^'^ M© yi~^ 
and Voc — 800 km s~^. 



uncertainty of ± 2000 K. Our modelling indicates that V 
sim = 40 km s^^ seems to better reproduce the obser- 
vation, especially the optical spectra. The gravity deter- 
mined by Villamariz et al. (,2002,1 gives a good fit of the 
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Fig. 6. Best fit (red dashed line) of the UV spectrum 
(black solid line) of HD 34078. The wind parameters 
are: M = IQ-^-^ Mq yr-\ v^o = 800 km 8"^ and 
logix/iboi = -7.0 
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Fig. 7. Zoom on the Iron line forests from Fig. |H| show- 
ing the good agreement between the predicted spectrum 
(dotted line) and the observed spectrum (solid line) and 
confirming the T^s estimate. 



Balmer lines, so that we adopt logg — 4.05 The other 
stellar parameters are gathered in Table El 

Figs. El and El show the fit of the Ha line and UV spec- 
trum of HD 34078. As we have shown that X-rays seem 
to be important for weak winds (see also next stars) and 
as HD 34078 shows no sign of a strong wind, we have 
adopted logLx/^boi = —7.0. Indeed, no X-ray measure- 
ment exists for HD 34078 and we have thus adopted the 
classical value for O stars (e.g. Chlebwoski & Garmany 
I1991[l . A reasonable agreement between the two types of 
mass loss indicators (Ha and UV lines) is found for M= 
10~9-^ M0 yr~^and a terminal velocity of 800 km s~^. 
Again, our value of M is lower than previous determina- 
tions (see Table OJ. 

Surprisingly, the derived terminal velocity is similar 
to or even lower than the escape velocity (1043 km s~^). 
However, given the large error on M and R, the escape 
velocity is also very uncertain. Also, a value of v^o lower 
than Vesc is possible since the escape velocity quoted here 
is the photospheric escape velocity, and a velocity in the 
wind of the order is obtained only in the outer at- 
mosphere where the local escape velocity is much lower. 
Moreover, the weakness of the wind features may actually 
lead to underestimate of the terminal velocity (see also 

sect. rm^ . 

For HD 34078, we have used the He and CNO abun- 
dances of Villamariz et al. H2002|l . They are nearly solar, 
except for C which is found to have an abundance of ^ 
1/2 solar. 



4.3. HD46202 

HD 46202 is an 09 V star situated in the Rosette nebula. 
An effective temperature of 33000 K gives the best fit of 
the optical He lines as shown in Fig. jSl As for HD34078, 
the too deep C ii AA6578,6582 doublet may indicate a 
slightly too low Teg: once again, increasing T^s by 1000 K 
improves the fit but gives a worse fit of the He lines, so 
that we think the value of 33000 K is reasonable within 
its uncertainty of ± 2000 K. A gravity of log g = 4.0 gives 
a good fit of the H7 fine (see Fig. jSJ . 

Figures El and El show the fit of the wind sensitive lines 
from which we derive a mass loss rate of 10~®'^ Mq yr~^ 
and a terminal velocity of 1200 km s~^. According to the 
X-ray detection, we have chosen log — —6.10 (see 
Table 01 . If X-rays are not included, a value of M as low 
as 10""'^° Mq yr~^ is required to fit the wind part of C IV 
AA1548,1551. Note that the core of the Ha line is stronger 
in the model, but as the observed profile seems to be some- 
what contaminated (possibly by a small nebular contribu- 
tion), we did not try to fit this core. As C iv AA1548,1551 
is the main M indicator and as log is quite high for 
this star, we have run test models including the high ion- 
isation states C V and C VI to check if the C ionisation 
was modified. They show that the C ionisation is indeed 
slightly increased, which implies to increase M by a factor 
of ~ 2 in order to fit C iv AA1548,1551. Hence, given the 
uncertainty in log (due to both uncertainties in Lx 
and iboi), we think this effect is negligible compared to 
other sources of errors for the M determination (see Sect. 
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Fig. 8. Best fit (red dashed line) of the optical spectrum 
(black sohd hne) of HD 46202. Here, T^tf = 33000 K, logy 
= 4.0, V sini = 30 km s"^ M = lO"^-^ Mq yr^^ and 
Voo — 1200 km s^^. The observed core of Ha is likely 
contaminated by small interstellar emission. 
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Fig. 9. Best fit (red dashed line) of the UV spectrum 
(black solid line) of HD 46202. The wind parameters are: 
M = IQ-^-^ Mq yr-i and v^o = 1200 km s~^ X-rays are 
included so that logLx/L^oi = —6.05. The lUE spectrum 
below 1200 A is not shown since the low S/N ratio does 
not allow any reliable comparison. 



Fig. 10. Best fit (red dashed line) of the observed He and 
Ha fines (black solid line) of HD 93028. The effective tem- 
perature is 34000 K, logg = 4.0 and V smi = 50 km s~^. 

E)) . We have also run test models for which the X-ray tem- 
perature was increased from 3 10^ K to 7 10^ K. Fitting 
C IV AA1548,1551 with this new X-ray temperature re- 
quired a slight increase 0.3 dex) of the mass loss rate. 
All previous studies give higher values of M (see Table OJ. 
As for Voo, the range of values derived by other authors 
is quite large and encompasses our estimate. This shows 
the difficulty of deriving v^o from weak wind line profiles. 
The low terminal velocity will be discussed in Sect. 17.211 

4.4. HD93028 

HD 93028 has a spectral type 09V and belongs to the 
young cluster CoUinder 228 in the Carina nebula. A value 
of V s'mi of 50 km s~^ was deduced from optical lines 
fits (Fig. and from previous studies (Table |21). The 
effective temperature we derive from He optical lines is 
34000 K. 

From the C iv AA1548,1551 line, we derive a termi- 
nal velocity of 1300 km s^^, slightly lower than other 
estimates (see Table O. We find that a mass loss rate 
of 10-9-5 M© yr^i gives a good fit of the far UV, UV 
and Ha spectrum without X rays. However, as we have 
shown previously. X-rays influences strongly the determi- 
nation of M in stars with weak winds (see Sect. 14. ll 14. 3() . 
Hence, although there is no measurement of X-rays for 
HD 93028, we adopted the classical value log = —7.0 
(Chlebowski & Garmany 11991(1 and then derived M — 
10^^ ° Mq yr~^ as shown in Figs. [TUI and [TTI The core of 
Ha is a little too strong in our best fit model, but the ob- 
served line shows evidences of interstellar contamination, 
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Fig. 11. Best fit (red dashed line) of the UV spectrum 
(black solid hue) of HD 93028. For this model, M = lO"^ " 
M0 yr~^, Woo = 1300 km s~^and we adopted log = 
-7.0. 

which is natural in a star forming region (see also the Ha 
profile of HD 93146). The only previous determination of 
mass loss rate for HD 93028 was made by Howarth & 
Prinja ifTM^ who found M = 10""^ M© yr ^, more than 
two orders of magnitude higher than our value. 
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Fig. 12. Best fit of the optical spectrum (black solid line) 
of HD 152590. The effective temperature is 36000 K, log 5 
= 4.1 and V sini — 66 km s~^. The blue dotted line is for 
a standard model {vturb — 20 km s^^ in the computation 
of the atmospheric structure and vturt from 5 to 175 km 
s~^ for the spectrum) while the red dashed line is for a 
model with vturb ~ 10 km s^^ (atmospheric structure) 
vtMrb = 10-175 km s~^ for the spectrum and additional 
metals (Ne, Ar, Ca and Ni). We see that this improved 
model leads to better fits of the He I singlet lines, leaving 
all other lines basically unchanged. 



4.5. HD152590 

HD 152590 is an 07.5Vz star. The distance estimate is 
difficult since it membership to Trumpler 24, Sco OBI or 
NGC 6231 is not completely established. Given the un- 
certainty in the distance, we simply adopt the mean value 
(see Table CJ. 

We adopted V smi = 66 km s^^ from Penny 1)1996(1 . 
The optical spectrum shown in Fig.^lis correctly repro- 
duced with an effective temperature of 36000 K. Note that 
initially, we had a problem to reproduce the He I singlet 
lines which were too weak in our models wheras all other 
lines were very well reproduced. This problem has been re- 
cently noted by Puis et al. (2005) when they put forward a 
discrepancy between CMFGEN and FASTWIND for these 
lines between 36000 and 41000 K for dwarfs. However, a 
more complete treatment of line blanketing appeared to 
solve this problem. Indeed, if we reduce the microturbu- 
lent velocity from 20 to 10 km s^^ in the computation 
of the atmospheric structure AND if we add some more 
species (Neon, Argon, Calcium and Nickel) we greatly im- 
prove the fit of the He I singlet lines without modifying 
the strength of other H and He lines (Hillicr et al. 120031 
already noted that the He I singlet lines were much more 



sensitive to details of the modelling than the triplet lines) . 
This is shown in Fig. Note that increasing the micro- 
turbulent velocity from 5 to 10 km s^^ in the compu- 
tation of the spectrum changes only marginally the line 
profiles. Hence, we attribute the origin of the discrepancy 
pinpointed by Puis et al. H2005|l to a subtle line-blanketing 
effect in this particular temperature range, and concern- 
ing only the He I singlet lines ^. Note that reducing vturb 
without including additional metals strengthens the sin- 
glet lines, but not enough to fit the observed spectrum. 
Hence the additional line-blanketing effects of Ne, Ar, Ca 
and Ni, although small (most lines are unchanged) is cru- 
cial to fit the He I lines around Tcff = 36000 K. Note that 
we usually restrict ourselves to models with vturt = 20 km 
and no Ne, Ca, Ar or Ni since the computational time 
is much more reasonable. For HD 152590, we found that a 

^ The problem occurs only for the He I singlet lines which 
have the 2p lp° state as their lower level. Thus the singlet 
problem is most likely related to the treatment of blanketing 
in the neighbourhood of the He I resonance transition at 584 
A. Detailed testing by F. Najarro & J. Puis (private commu- 
nication) supports these ideas. 
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Fig. 13. Fits of the observed Ha line (black solid line) of 
HD 152590 for a model with M = IQ-'^-'^^ Mq jt'^ (red 
dashed line) and for a model with M — 10~®'^^ Mq yr~^ 
(blue dotted line). The terminal velocity in both models 
is 1750 km s~^. Note the insensitivity of Ha line profile to 
the mass loss rate. 

gravity log 5 = 4.10 ± 0.1 gives the best fit of the Balmcr 
lines (in particular H7) 

The terminal velocity of HD 152590 estimated from 
C IV AA1548,1551 is 1750 km s^^, in good agreement or 
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Fig. 14. Same as Fig. II 31 for the UV range. 



Fig. 15. Best fit (red dashed line) of the observed He and 
Ha lines (black solid line) of HD 93146. The effective tem- 
perature is 37000 K, logg = 4.0 and V smi = 80 km s~^. 



lower than previous estimates (Table Ojl- The estimate of 
the mass loss rate is much more difficult for this star. In 
fact, we have not been able to fit simultaneously the UV 
lines and Ha. If the former are correctly reproduced (with 
M = 10~^'^^ M© yr~^), then the later has a too strong ab- 
sorption in its core, and if Ha is fitted (with M = 10~^'^^ 
M0 yr~^), C IV AA1548,1551 is too strong. This is shown 
in Fig. El and We have tried without success to in- 
crease the (3 parameter to improve the fit (an increase of 
(3 leading to a weaker Ha absorption) . The fits of Fig. [T^ 
and 1141 are for /3 = 1.2 and even for this quite high value 
for a dwarf star, the Ha core is not perfectly reproduced. 
A possible explanation is the presence of a companion for 
HD 152590 (Gieseking|1982). In that case. Ha may be di- 
luted by the continuum of this secondary whereas the UV 
spectrum may be unaffected provided the companion is a 
later type star than HD 152590 without strong UV lines. 
However, adopting a conservative approach, we adopt the 
Ha mass loss rate (lO"^'"^® M© yr~^) as typical, keeping 
in mind that it may well be only an upper limit. 

4.6. HD93146 

HD 93146 is an 06.5V((f)) star in the Carina nebula and 
belongs to the cluster Cr 228. 

We adopt V sini = 80 km s^^ from our fits and previ- 
ous determinations (see Table 12) • Fig. 1151 shows our best 
fit to the He optical spectrum between 5000 and 6000 A 
for which an effective temperature of 37000 K is derived. 
Notice that this fit is not perfect, but it is actually the best 
we could get. Increasing Toff may help reduce the He I ab- 
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Fig. 16. Best fit (red dashed line) of the UV spectrum 
(black sohd line) of HD 93146. For this model, M = 
10-7-25 Mq yr-iand Voc = 2800 km s^i. 
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Fig. 17. Best fit (red dashed line) of the observed He and 
Ha fines (black solid line) of HD 42088. The effective tem- 
perature is 38000 K, logg = 4.0 and V smi = 60 km s"^. 



sorption, but it increases too much the He II strength. 
Moreover the UV photospheric lines are very well repro- 
duced with this Tcs (see Fig. I16|l . As we do not have re- 
liable gravity estimators, we assume logg = 4.0 since this 
value is typical of dwarfs (Vacca et al. 119961 Martins et al. 

Fig. ^1 shows our best fit of the (far) UV spectrum 
of HD 93146. The terminal velocity is 2800 km s'^ and 
the mass loss rate is lO-^'^^ Mq yr'-^. For higher values, 
N IV A1718 displays a too strong blueshifted absorption. 
The Ha profile of Fig. El confirms partly this value of 
M since the line is correctly reproduced, under the uncer- 
tainty of the exact depth of the core which is contaminated 
by nebular emission. Previous estimates are in failry good 
agreement with the present one (Table E)). 

4.7. HD42088 

HD 42088 is a 06.5 V star associated with the H II region 
NGC 2175. It also belongs to the class of Vz stars. Note 
that the distance to this star is poorly constrained so that 
its luminosity is the least well known of all stars of our 
sample. The rotational velocity is chosen to be 60 km s^^ 
in view of the determinations of Penny H1996|l - 62 km s^^ 
- and Howarth et al. [1997) - 65 km s"'^. The fit of optical 
He lines above 5000 A leads to an estimate of the effective 
temperature which is found to be ~ 38000 K as shown 
by Fig. El This estimate also relies on the fit of UV lines 
since the number of optical indicators is small. We adopt 
log (7 = 4.0 (from Vacca et al. I1996|l since we do not have 
strong gravity indicators. 
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Fig. 18. Best fit (red dashed line) of the UV spectrum 
(black solid fine) of HD 42088. For this model, M ^ 10"* 
Mq yr-^and v^o = 1900 km s~^. 



The terminal velocity is derived from the blueward ex- 
tension of the absorption in C iv AA1548,1551 and is 1900 
km s"^. Previous determinations go from 2030 km s^^to 
2550 km s""'^Given the fact that we adopted a microtur- 
bulent velocity of 190 km s-^in the outer wind (10 % of 
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Fig. 20. Best fit of the UV spectrum (black solid line) of 
HD 93204. For this model, M = IQ-^ ^s j^q yr^^and v^o 
= 2900 km s"^. The red dashed line is a clumped model 
with foe =0.1 and the blue dotted line is a homogeneous 
model. 



Woo), the absorption actually extends up to 2100 km s^-'^in 
the model, in good agreement with other determinations. 
Concerning the mass loss rate, it turns out that a value 
of 10"^ M0 yr^^gives a reasonable fit of the main UV 
lines and Ha, although for the latter the very core is not 
correctly fitted but may suffer from nebular contamina- 
tion (see Fig. [T7|) . The best fit model is shown in Fig. ITSl 
Our mass loss rate determination based on both Ha and 
UV lines gives a much lower value than ever found for 
this star (Table OJ). But the UV lines produced by models 
with mass loss rates much higher than our adopted value 
are much too strong compared to the observed spectrum, 
forcing us to adopt such a low M . 

4.8. HD93204 

HD 93204 (05V((f))) is a member of the young cluster 
Trumpler 16 in the Carina complex. We adopt the value 
130 km s^^ for V sini in our fits, which helps to derive an 
effective temperature of 40000 K (see Fig. ll9f) A gravity of 
log 5 = 4.0 is compatible with the observed Balmer lines. 

Fig. 1201 shows the fit of the UV spectrum. To fit rea- 
sonably all the UV lines, we had to use clumped models. 
This is especially true for O v A1371 since as previously 
shown by Bouret et al. H2003I I2005|l this line is predicted 
too strong in homogeneous models. In our case, the use 
of clumping with the law given in Sect. 13.21 and foo = 0.1 
improves the fit of O v A1371 as weU as N iv A1718, as 
shown in Fig. [201 Reducing Toff does not solve the prob- 
lem since in that case O V A1371 is weaker but N iv A1718 



gets stronger. We derive a mass loss rate of lO^^'^^ Mq 
yr^^ and a terminal velocity of 2900 km s^^. Due to the 
high level of nebular contamination of Ha, we can not use 
this line to constrain M (see Fig. I19f) . Our value of M is 
slightly smaller (factor 4) than that of Howarth & Prinja 
(1989) (Table Ell mainly due to the inclusion of clumping 
in our models. As for Voo, our estimate is well within the 
range of values previously derived (see Table EJ. 



4.9. HD15629 

HD 15629 is classified as 05V((f)) and belongs to the 
star cluster IC 1805. The projected rotational velocity is 
found to be 90 km s^^ by several authors and we adopted 
this value which give a good fit of optical and UV pho- 
tospheric lines. The optical spectrum presented in Fig. 1211 
indicates an effective temperature of 41000 K. This is in 
good agreement with the recent determination of Repolust 
et al. iPTOjl who found 40500 K. We adopted log 5 = 3.75 
since it gives a reasonable fit of Balmer lines fFig.l21|l and 
it is close to the value derived by Repolust et al. H2004|l 
who derived logg = 3.70. 

The best fit model of the UV spectrum is shown in 
Fig. 1221 and Ha is displayed in Fig. |221 The main param- 
eters for this model are M = lO'^-^ Mg yr'^ Voa = 2800 
km s"^ and foo = 0.1. We also show on these figures a 
model without clumping and with the mass loss rate of 
Repolust et al. which is higher - IQ-^-^^ Mq yr'^- 

than our derived value. Once again the inclusion of clump- 
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Fig. 21. Best fit (red dashed line) of the observed He and Fig. 22. Best fit (red dashed line) of the observed Ha line 
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H hnes (black solid hne) of HD 15629. The effective tern- (black sohd line) of HD 15629. For this model, M = 10 
perature is 41000 K, logg — 3.75 and V smi = 90 km yr~^, Voo — 2800 km s^^ and /oo — 0.1. We also show 

. a model with M — io~^-^^ M© yr"^ and no clumping 

(blue dotted line) as derived by Repolust et al. ((20031 • 



ing is necessary to correctly reproduce both O v A1371 
and N IV A1718. With the Repolust et al. pTOIjl M and 
no clumping, CNO abundances have to be reduced by a 
factor of 3 to give reasonable fits, and even in that case 
the O V A1371 line is too strong. Such a reduction of the 
abundances is unlikely for a Galactic star. For our best 
fit, we have adopted the CNO solar abundances recently 
claimed by Asplund H2004|l since they are slightly lower 
than those of Grevesse & Sauval 1^1998 ) and allow a fit of 
the UV lines with a slightly higher (0.25 dex) mass loss 
rate compared to the later values. Note that in our final 
best fit, the core of Ha is not perfectly fitted. However, 
we suspect that the strange squared shape of the observed 
line core is probably contaminated by weak nebular emis- 
sion. In support of the nebular contamination we note the 
following: if we adopt the mass loss rate of Repolust et 
al. 1)2004(1 . the fiux level in the line core is correct, but 
the line is slightly narrower in the remainder of the profile 
compared to the observed profile, while with our Af , the 
line is well fitted except in the very core. Increasing the 
flux level in the core in models with our M requires the 
adoption of fi — 1.7 which is high for a dwarf. In that 
case again, although the flux level in the core is correct, 
the synthetic line profile is too narrow. We are then rather 
confident that the observed line core is somewhat contam- 
inated and that our mass loss rate is correct. The use of 
clumping explains partly the discrepancy with the result 
of Repolust et al. 1(2004(1 . Concerning the terminal veloc- 
ities, previous extimates range from 2810 to 3220 km s^^ 
in reasonable agreement with our value. 




Mi) 

Fig. 23. Best fit (red dashed line) of the UV spectrum 
(black solid line) of HD 15629. The wind parameters are 
the same as in Fig. 122 The lUE spectrum below 1200 A 
is not shown since the low S/N ratio does not allow any 
reliable comparison. The blue dotted line is the model with 
the Repolust et al. 1(2004(1 parameters: note the too strong 
O V A1371 and N iv A1718 lines. 
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Fig. 24. Best fit (red dashed line) of the observed He and 
Ha lines (black solid line) of HD 46223. The effective tem- 
perature is 41500 K, logg = 4.0 and V sini = 130 km s~^. 

4.10. HD46223 

HD 46223 belongs to the Rosette cluster (NGC 2244) and 
has a spectral type 04V((f+)). A projected rotational ve- 
locity of 130 km s^^ was adopted from the fit of optical 
lines. The upper panels of Fig. [31 show the fit of He op- 
tical lines with a model for which Tcs — 41500 K. Note 
that this effective temperature also gives a reasonable fit 
of the UV spectrum (see Fig. I25|l . The subsequently de- 
rived stellar parameters are gathered in Table El As we do 
not have reliable gravity indicators, we adopt logg = 4.0 
As regards the terminal velocity, we find Voo = 2800 
km from the UV resonance lines. This is in fairly 
good agreement with previous estimates (see TablejSJ- The 
mass loss rate is derived from Ha and the UV resonance 
hues. The adopted value for M is 10"*^-^ Mq yr~^ As 
for HD 93204, clumping was necessary to fit O v A1371 
and N iv A1718. Since the inclusion of clumping leads to 
mass loss rates lower than in homogeneous winds, this ex- 
plains partly why our estimate is nearly a factor 5 lower 
than most previous estimates for this star which did not 
use clumping (see Table O- Note that in our models, the 
inclusion of clumping reduces the strength of N v A 1240 
which is then less well fitted than in the case of the ho- 
mogeneous model. However, the very blue part of the ab- 
sorption profile is contaminated by interstellar Lyman ab- 
sorption rendering the exact line profile uncertain. 

4.11. HD93250 

HD 93250 is a well studied O dwarf of the Trumpler 16 
cluster in the Carina region. It is a prototype of the re- 
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Fig. 25. Best fit of the UV spectrum (black solid line) of 
HD 46223. For this model, M = 10"^-^ Mq yr^^and Voo 
= 2800 km s"'^. The red dashed line is a clumped model 
with foe =0.1 and the blue dotted line is a homogeneous 
model. 



cently introduced 03.5 subclass (ST 03.5((f+)) Walborn 
et al. 2002). 

Optical lines indicate a projected rotational velocity 
of 100 km and an effective temperature of ^ 46000 
K (mainly from the strength of He i A4471, see Fig. 1^ . 
However, Fe line forests in the UV are more consistent 
with a value of 42000-44000 K as displayed in Fig. OHl 
For such a Toff He i A4471 is a little too strong in the 
model. However, this seems to be the case of all H and He 
optical lines, possibly due to the fact that HD 93250 may 
be a binary (see Repolust et al. I2004|l which may also 
be advocated from the fact that the absorption of C iv 
AA1548,1551 is not black despite the strength of the line 
(allowing the study of discrete absorption components). 
Hence, we rely mainly on the UV and we adopt a value of 
44000 K for the effective temperature of HD 93250. This 
value is in reasonable agreement with the determination 
of Repolust et al. who found 46000 K. We adopted 

logg = 4.0 from Repolust et al. (|2004p and our fit of H7. 
Note that the estimated mass for this star is especially 
high (144 Mq), which may make HD 93250 one of the 
highest mass stars known. However, the uncertainty on 
the mass determination is huge, and HD 93250 is also 
suspected to be a binary. Hence, we caution that the mass 
given in Table El is only indicative. 

The determination of the wind parameters relies on 
Ha and on several strong UV hues: N v A1240, O iv 
AA1339,1343, O v A1371, C iv AA1548,1551, He 11 A 1640 
and N iv A1718. The terminal velocity deduced mainly 
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I1989f) . However, we use a microturbulent velocity of 200 
km in the outer part of our model atmosphere for this 
star, so that in practice, the absorption extends up to 3200 
km s~^. As regards the mass loss rate, we actually found 
that it was impossible to find a value for M which would 
produce reasonable fits of all UV lines in homogeneous 
winds. Indeed, O v A1371 was always too strong and N iv 
A1718 too weak. Reducing the effective temperature does 
not improve the situation, since values as low as 40000 K 
are required to fit O V A1371, and in that case the other 
UV lines are not correctly fitted so that again, we had to 
include clumping. In the end, we find that a mass loss rate 
of 10^^-^^ M© yr^-^with a clumping factor foo = 0.01 gives 
a reasonable fit, as displayed in Fig. |23 This value of M 
is lower than the determination of Repolust et al. I)2004|l 
- 10~^'^^ Mq yr~^- relying only on Ha. We will return 
to this in Sect. 17.2'!^ Note that the value of foo we derive 
is quite small, but not completely unrealistic in view of 
recent results presented by Bouret et al. (^2005^) indicating 
foo = 0.02 and 0.04 for two 04 stars. 



Fig. 26. Best fit (red dashed line) of the observed He and 
Ha lines (black solid line) of HD 93250. The effective tem- 
perature is 46000 K, log5 = 4.0, V sini = 110 km s'^M 
= 10-6-25 yr-i^ = 3000 km s^^ and foe = 0.01. 



from C IV AA1548,1551 is 3000 km s^^, shghtly lower than 
the previously derived values which are between 3250 km 
(Repolust et al. r2004|l and 3470 km s-^ (Bernabeu 
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Fig. 27. Best fit of the UV spectrum (black solid line) of 
HD 93250. For this model (red dashed hue) M = IQ-^ ^s 
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3000 km s-^and foo = 0.01. 
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Fig. 28. Determination of effective temperature from UV 
Fe line forests. Solid line is the observed spectrum, dotted 
line a model with Toff ~ 40000 K, long dashed line a model 
with Tcff — 42000 K, short dashed line a model with TeS 
= 44000 K and dot-dashed hue a model with Tcff = 46000 
K. See text for discussion 



5. Role of X-rays and magnetic field in weak-wind 
stars 

Several of our sample stars have published X-ray fluxes. 
Chlebowski & Garmany (T9!JI) report X-ray measure- 
ments for HD 38666, HD 46202, HD 152590, HD 42088 
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and HD 46223, while Evans et al. ^M^i give X-ray lu- 
minosities for HD 93204 and HD 93250. These high en- 
ergy fluxes may have important consequences on the at- 
mosphere structure since, as shown by MacFarlane et al. 
P994|l . the ionisation fractions may be significantly al- 
tered. These authors also demonstrated that the effect of 
X-rays was higher in low-density winds: ionisation in early 
O stars is almost unchanged by X-rays, while in early B- 
stars changes as large a factor 10 can be observed between 
models with and without X-rays. The reason for such a be- 
haviour is that 1) X-rays produce higher ionisation state 
through single ionisation by high energy photons and the 
Auger process and 2) the ratio of photospheric to X-ray 
flux decreases when effective temperature decreases, im- 
plying an increasing role of X-rays towards late type O and 
early B stars (see MacFarlane et al. 11994(1 . Moreover, the 
lower the density, the lower the recombinations to com- 
pensate for ionisations so that we expect qualitatively an 
even stronger influence of X-rays in stars with low mass 
loss rate. Since some of our sample stars are late type O 
stars with low density winds. X-rays can not be discarded 
in their analysis. Indeed, the Carbon ionisation fraction - 
and thus the strength of the C iv AA1548,1551 line and 
the derived mass loss rates - can be altered. 

In this context, we have first run test models for HD 
46202 and HD 93250. For HD 93250, the inclusion of X- 
rays did not lead to any significant change of the ionisation 
structure as expected from the above discussion. Indeed, 
the main wind line profiles were not modified (see also 
MacFarlane et al. HM^ Pauldrach et al. [T!M1 IMTTjl . in- 
dicating that X-rays are not crucial for the modelling of 
these lines in such high density winds. Of course, other 
lines are well known to be infiuenced by X-rays (e.g. O vi 
AA1032,1038) but are not used in this study to derive the 
stellar and wind parameters. However, in the case of HD 
46202 the ionisation structure in the wind is strongly mod- 
ified which leads to a weaker C iv AA1548,1551 line (for 
a given M) as displayed in Fig. I^Hl Indeed, the ionisation 
fraction of C IV is reduced: this is displayed in Fig. 1201 
where we see that in the model giving the best fit with- 
out X-rays, C iv is the dominant ionisation state, while 
when X-rays are included, it is no longer the case. Fitting 
C IV AA1548,1551 thus requires a higher mass loss rate. 
In practice, the change in the C iv AA1548,1551 profile 
when X-rays are included is equivalent to a reduction of 
the mass loss rate by a factor of ~ 10 in models without 
X-rays. Given this result, we have included X-rays in our 
modeUing of the atmosphere of HD 38666, HD 46202, HD 
34078 and HD93028. For the two former stars, X fiuxes 
from the literature have been used while for the two latter 
ones, we simply adopted log = —7.0. 

A question which remains to be answered concerning 
the X properties of such weak wind stars is the origin 
of the X-ray emission. Indeed, it is usually believed that 
shocks in the wind due to instabilities in the line driving 
mechanism are responsible for the production of such high 
energy photons (Lucy & White 119801 Owocki, Castor & 
Rvbicki . 1988,1 . This scenario seems to apply to the strong 
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Fig. 29. Effect on X-rays on the C iv AA1548,1551 line. 
The observed profile is the solid line, the initial model is 
the dotted line and the model with X-rays and the same 
M is the dashed line. See text for discussion. 



wind star Zeta Pup (Kramer et al. 12003(1 . However, recent 
observations by Chandra have revealed that for the BOV 
star T Sco and for the Trapezium stars, most of the lines 
emitted in the X-ray range were too narrow to have been 
produced in the wind up to velocities of the order V(yQ cis ex- 
pected in the wind-shock scenario (see Cohen et al. IM)^ 
Schulz et al. 12003). And these lines are also not formed 
very close to the photosphere as predicted by a model 
in which the X-ray emission is due to a hot corona (e.g. 
Cassinelli & Olson [T979(l . Actually, such lines are more 
likely to be formed in an intermediate region. This may be 
explained in the context of magnetically confined winds: 
in this scenario, the presence of a magnetic field confines 
the outflow and channels it into the equatorial plane where 
shocks produce X-ray emission above the photosphere but 
not in the upper atmosphere (See Babel & Montmerle 
11997(1 . This model has been recently refined by Ud'Doula 
|2002(l who have investigated the structure of 



& Owocki 

both the wind outfiow and the magnetic field through time 
dependent hydrodynamic simulations. In particular, they 
estimated from simple arguments the strength of the mag- 
netic field required to confine the wind (hereafter B') and 
thus to lead to shocks in the equatorial plane. 

In Table01 we have gathered different properties of the 
stars of our sample showing X-ray emission: the X-ray lu- 
minosity (ix), the mechanical wind luminosity (Lwind = 
iMu^) for our M and M from Vink et al. ^UU^, and 
B' again considering our derived M and Vink's M. We 



see that for weak winds, 



becomes of the order unity 



which shows the increasing importance of X-rays as the 
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wind becomes less and less dense. In addition, Table ^ 
shows that the magnetic field strength required to confine 
the wind is low for weak-wind stars, showing the increas- 
ing role of magnetic field when M decreases. Given these 
results and the above discussion, we may speculate that 
our weak wind stars may have magnetically confined winds 
(although no detections of magnetic field exist for them). 
In that case, one may wonder how our results would be 
modified. Fig. 8 of Ud'Doula & Owocki pIT)^ shows that 
the mass flux {pv) is reduced close to the pole and en- 
hanced near the equator, but their Table 1 reveals that 
the total mass loss is only reduced by a factor < 2 even in 
the case of strong confinement. Hence, using classical ID 
atmosphere models should lead to correct values for the 
mass loss rates within a factor of two, even if magnetic 
confinement exists. 

As regards this last point, a comment on the shape 
of the line profiles is necessary. Indeed, the only O stars 
with a detected magnetic field - Thetal Ori C (Donati et 
al. 12002 ) - shows unusual features which may be related to 
the geometry of magnetically confined wind (Walborn et 
al. 119851 Gagne et al. . 1200511 . The absence of such unusual 
features in the spectra of our stars with weak winds may 
argue against such a confinement. However, as Thetal Ori 
C is the only example and as there is no theoretical pre- 
diction of the change in the shapes of wind lines in the 
presence of magnetic confinement - which also probably 
depends on the tilt angle between the nagnetic and rota- 
tion axis -, we can not completely rule out the existence 
of magnetic confinement in our weak wind stars. 

6. Sources of uncertainty for the M determination 

In this section, we investigate the various sources of un- 
certainty of our determinations of mass loss rates both on 
the observational side and on the modelling side. 

6.1. Observational uncertainties 

Under the term "observational uncertainty" , we gather all 
the efi:ects which can influence the shape of the observed 
line profiles, especially Ha. The first source of uncertainty 
is the S/N ratio. However, in most of the stars studied 
here, this ratio is good (~ 100) and does not affect the 
analysis. The second source of uncertainty comes from the 
normalisation of the spectra. This is a general and well 
known problem which can affect the strength of lines, es- 
pecially in the case of weak lines. In our spectra the main 
difficulties arise in the N v A1240 and Ha regions. For the 
former, this is due to the presence of the broad Lyman 
a absorption around 1216 A which renders uncertain the 
exact position of the continuum. We simply check that the 
strength of the emission part of the profile in the models 
is on average consistent with the observed line, leaving 
aside the bluest part of the absorption. The case of Ha 
is more critical. The normalisation can be hampered by 
the S/N ratio: a low ratio will not allow a good identifica- 
tion of the continuum position. The use of echelle spectra 
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Fig. 30. Effect of the inclusion of X-rays on the Carbon 
ionisation structure. Shown is the ratio of the population 
of C ionisation state i to the total H population as a func- 
tion of Rosscland optical depth. The red dashed line corre- 
sponds to the model giving the best fit to the UV spectrum 
when X-rays are not included, while the black solid line 
is the best fitting model with X-rays. In the former case, 
C IV is the dominant ionisation state in the wind, while 
it is no longer the case in the latter model (C V being 
the dominant). Note that Cothcr refers to the difference 
between the total C population and the sum of C ii, C ill 
and C IV, and is mainly C v here. 



renders also difficult the identification of the continuum 
since the wavelength range around the line of interest in a 
given order is limited to ~ 60 A. We estimate that taken 
together, these effects induce an uncertainty < 0.02 on 
the absolute position of the Ha core. Of course. Ha is also 
contaminated by nebular emission. When present, such an 
emission precludes any fit of the very core of the line. But 
the high resolution of our spectra allows a fit of ~ 80—90% 
of the stellar profile, excluding the very core. 

6.2. Photospheric Ha profile 

Our estimates of M rely on the fit of both the UV wind 
sensitive lines and Ha. In low density winds, Ha is essen- 
tially an absorption profile for which only the central core 
is sensitive to mass loss rate. In order to derive reliable 
values of M it is thus important to know how robust the 
prediction of the photospheric profile is, since it will dom- 
inate over the wind emission. This is of much less impor- 
tance in high density winds where the lines are dominated 
by wind emission. 
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Table 4. X-ray properties of our sample stars with known X-rays fluxes. Lx is from Chlebowski & Garmany H1991|l 
for HD 38666, HD 46202, HD 152590, HD 42088 and HD 46223, and from Evans et al. ipnUSjl for HD 93204 and HD 
93250. iwind is the mechanical wind luminosity. Values with "Vink" are those for which M is taken from Vink et al. 
H2001|> mass loss recipe. B' is the value of the magnetic field for which confinement begins (corresponding to jy^, = 1 
in the formalism of Ud'Doula & Owocki O002|l . 



HD 


log Lx 


log-^ 


log Lwind 


1„„ T Vink 
iOg ^wind 




loir 
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B' 






[erg s"^] 




[erg s"^] 


[erg s"^] 






[G] 


[G] 


38666 


31.37 


-6.87 


32.16 


34.25 


-0.79 


-2.88 


7 


75 


46202 


32.40 


-6.05 


32.76 


34.43 


-0.36 


-2.03 


11 


72 


152590 


32.51 


-5.86 


34.20 


34.83 


-1.69 


-2.32 


60 


125 


42088 


32.38 


-6.43 


34.06 


35.89 


-1.68 


-3.51 


33 


269 


93204 


32.06 


-7.03 


35.67 


36.31 


-3.61 


-4.25 


137 


286 


46223 


32.62 


-6.53 


35.89 


36.42 


-3.27 


-3.80 


180 


332 


93250 


33.22 


-6.53 


36.20 


37.20 


-2.98 


-3.98 


153 


485 



To check the CMFGEN prediction in a low density 
wind, we have compared the Ha line with that predicted 
by FAST WIND, the other non-LTE atmosphere code in- 
cluding wind and line-blanketing widely used for optical 
spectroscopic analysis of massive stars (see Santolaya-Rey 
et al.CP?, Puis et al. lMlH| . The test model was chosen 
with the following parameters: TcS — 35000 K, logg — 4.0, 
Af = 10"^ M© yr"^ and /3 = 0.8. This set of parameters 
is typical of the stars with weak winds analysed in the 
present study, and Ha should not be too much contam- 
inated by wind emission. The result of the comparison 
between CMFGEN and FASTWIND is given in Fig. EH 
We see that the agreement between both codes is very 
good. This is not a proof that the predicted profile is the 
correct one, but it is at least a kind of consistency check. 

We have also investigated another effect which can al- 
ter the shape of the Ha line core: the number of depth 
points included in the models. Indeed, the thinner the spa- 
tial sampling, the better the line profile. This means that 
a too coarse spatial grid should introduce errors in the 
determination of M from Ha. We have run a test model 
taking the best fit model for HD 34078 and increasing the 
number of depth points from 72 to 90: a thinner spatial 
grid leads to a slightly less deep line core, but the differ- 
ence is only of 0.01 in terms of normalised flux. This is 
lower than any other observational uncertainty (see Sect. 
I6.1|l so that we have adopted ~ 70 depth points in all our 
computations ^. 

In conclusion, there is no evidence that the photo- 
spheric Ha profile is not correctly predicted by our models. 

6.3. lonisation fraction 

In the low density winds, which correspond to late O type 
dwarfs in the present study, the final word concerning the 
mass loss rate is often given by C iv AA1548,1551. Indeed, 
Ha becomes almost insensitive to M in these cases, and 

^ choosing 90 depth points significantly increases the re- 
sources required for the computation. 




Ml) 

Fig. 31. Comparison between CMFGEN (red dashed fine) 
and FASTWIND (black solid line) Ha profile. The model 
is for Teff = 35000 K, logg = 4.0, M = IQ-^ M© yr"! and 
(3 = 0.8. The agreement between both codes is good. 



the other main wind sensitive UV line, N v A1240, is al- 
most absent from the spectra due to the reduced effective 
temperature. Other indicators such as Si IV AA1394,1403 
or N IV A1718 are still present, but they are weaker than 
C IV AA1548,1551 and become rapidly insensitive to any 
change of the mass loss rate. For more standard winds, al- 
most all indicators can be used together to derive M. We 
show in Fig.|23|the variation of the C iv AA1548,1551 line 
profile when the mass loss rate is decreased from 10~*'^ 
down to 10-^-^ Mq yr-i for the case of star HD 46202. 
We clearly see that C iv AA1548,1551 is still sensitive to 
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Fig. 33. Determination of M in low density winds (HD 46202). (a) shows the variation of the C IV AA1548,1551 hne 
profile when M is reduced (solid line: observed line; dash-dotted line: M = 10~^'^ Mq yr~^; long dashed line: M = 
10~®'^ Mq yi:~^; short-dashed long-dashed line: M = 10~^'^ Mq yr~^). Notice that the weaker the line, the more 
prominent the presence of photospheric lines (mostly from iron) superimposed to the C iv absorption, (b) shows the 
behaviour of Ha under the same changes. See text for discussion. 
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Fig. 32. Same as Fig. |2]but for H and He optical lines. 
The agreement between both codes is good. 



changes in M even for such low values. In parallel, we see 
that Ha is essentially unchanged in this regime of M. 



However, relying on only one line to assign final mass 
loss rates may be risky. We have highlighted in paper I that 
erroneous mass loss rates may be derived in the case where 
the C IV ionisation fraction is incorrectly predicted. This 
is still true here, since fitting the observed profile gives 
the right value of M x qc iv (stc iv being the ionisation 
fraction of C IV) but not necessarily the right M. In Fig. 
1341 we compare the ionisation fractions predicted by the 
CMFGEN best fit models to the values derived by Lamers 
et al. 1999] (hereafter L99) for dwarfs. The ionisations frac- 
tions are defined by 



qc IV 



_ /o^.2"'C iv{x)dx 
Io.2 ncix)dx 



(7) 



where x = and nc iv and nc are the number densities 
of C IV andtj respectively. At first glance, the CMFGEN 
ionisation fractions seem to be ~ 2 orders of magnitude 
higher than the L99 results, and that in spite of the few 
lower limits in the latter data. However, several comments 
can be made: 

- First, the work of L99 is based on previous mass loss 
rate determinations, mainly from Ha (Puis et al. I1996L 
Lamers & Leitherer I1993f) or from predictions for their 
dwarf subsample (Lamers & Cassinelli ll996|l . In the latter 
case, M is derived from the modified wind momentum - 
luminosity relation, so that any error in the calibration can 
lead to incorrect mass loss rate. Moreover, the uncertainty 
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of such a method due to the fact that a given star can 
deviate from a mean relation may introduce a bias in the 
derived ionisation fraction. Concerning the mass loss rates 
derived from Ha, Lamers & Leitherer ()1993|l use the line 
emission strength to determine M . However, in most O 
dwarfs Ha is in absorption so that the determination of 
the emission part of the line filling the photospheric profile 
may be uncertain. Puis et al. 1)1996(1 also use Ha to derive 
M but give only upper limits in the cases of thin winds. 
As L99 adopt these upper limits as the real values, we 
should expect the derived ionisation fractions to be lower 
limits. 



- Second, there is a significant shift in terms of pa- 
rameter space sampled by our results and that of L99: we 
have stars with 33000 < Toff < 44000 K and -17.3 < 
log < p X -14.4 (where < p >= j^^^ and ro.5 is 
the radius at which the velocity reaches half the termi- 
nal velocity) while L99 have 38000 < Toff < 50500 K and 
— 15 < —log < p >< —13.4, although both studies have 
stars of late and early O spectral types. Concerning effec- 
tive temperatures, part of the discrepancy comes from the 
use of line-blanketing in our models, which is known to 
reduce Tgs compared to unblanketed studies. But for den- 
sities, the explanation may again come from the fact that 
the adopted mass loss rates (and consequently the den- 
sities) in one or the other study are not correct. Can we 
discriminate between them? An interesting point is that 
3 stars are common to our study and that of L99: they 
are shown linked by dotted lines in Fig. |^ If we consider 
the fact that line-blanketing may explain the lower Tgff in 
our study, and the fact that for these stars the ionisation 
fractions derived by L99 are only lower limits, then the 
ionisation fractions predicted by CMFGEN are not neces- 
sarily too high. And if in addition we argue that our study 
investigates a density range not explored by L99, then we 
can not conclude that the ionisation fractions predicted 
by CMFGEN are wrong since no comparison can be made 
for very low mean densities. 



How could we test more strongly the wind ionisation 
fractions of our atmosphere models? One possibility is 
offered by the analysis of far UV spectra. Indeed, this 
wavelength range contains a number of lines formed 
in the wind from different ions of the same elements. 
Such a test will be done in a subsequent paper, based 
on FUSE observations of Vz stars in the LMC. But we 
can already mention that several studies of supergiants 
in the Magellanic Clouds using FUSE + optical data 
do not reveal any problem with the CMFGEN wind 
ionisation fractions, except that clumping must be used 
to reproduce a couple of (but not all) hues (see Crowther 
et alEOnS Hillier et al. IMH Evans et nl HUU^ . Hence 
in the following, we assume that the ionisation fractions 
given by CMFGEN are correct. 
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Fig. 34. C IV ionisation fractions in CMFGEN models 
(filled symbols) and from Lamers et al. 119991 (open sym- 
bols) for dwarfs as a function of effective temperature (up- 
per panel) and mean wind density (lower panel) . The dot- 
ted lines link objects in common between this work and 
the study of Lamers et al. ()1999,) . See text for discussion. 

6.4. Abundances 

Although our mass loss determination relies on both Ha 
and UV lines, we usually give more weight to the UV diag- 
nostics since the absorption profile of Ha can be shaped by 
other parameters than M {/3, clumping). But the UV lines 
depend more strongly on abundances than Ha. Hence, we 
have to estimate the error we make on the M determina- 
tion from UV lines due to uncertain abundances. We have 
already seen in Sect. 14.91 that adopting the CNO abun- 
dances of Asplund (l|2004p instead of those of Grevesse & 
Sauval 1(1998(1 - which corresponds on average to a reduc- 
tion by a factor of ~ 3/4 - leads to an increase of M 
by 0.25 dex. We have also run test models for a low lu- 
minosity star (HD 46202). It turns out that reducing the 
CNO abundances by a factor 2 implies an increase of the 
mass loss rate of the order of 2-2.5 in order to fit C IV 
AA1548,1551 since this line is not saturated in low density 
winds and thus its strength is directly proportional to the 
number of absorbers.. 

How different from solar could the CNO abundances 
of our sample stars be ? Given the estimated distances, 
it turns out that all stars are within 3 kpc from the the 
sun. Determinations of abundances through spectroscopic 
analysis of B stars (Smartt et al. 120011 RoUeston et al. 
1^00.) reveal the following gradients: —0.07 dex/kpc for C, 
O and Si, and -0.06 to -0.09 dex/kpc for N. Similarly, 
Pilyugin et al. ((2003(1 derive an Oxygen abundance gradi- 
ent of —0.05 dex/kpc from studies of HII regions. Taken 
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together, these results indicate that on average we do not 
expect variation of CNO abundances by more than ± 0.25 
dex for our sample stars. This means that adopting a so- 
lar metallicity leads to an error of at most 0.3 dex on the 
mass loss rate determination. 

Given the above discussion, we estimate the error on 
M due to uncertainties in the CNO abundances to be of 
~ 0.3 dex. 

6.5. Advection / adiabatic cooling 

In low density winds, two processes may affect the 
ionisation structure: advection and adiabatic cooling. 
The former is rooted in the fact that for low densities, 
the timescale for recombinations becomes longer than the 
timescale for transport by advection. Thus the ionisation 
structure can be significantly changed. The latter process 
(adiabatic cooling) lowers the temperature in the outer 
part of the atmosphere where the heating processes 
(mainly photoionisations) are less and less efficient due 
to the low density, implying also a modification of the 
ionisation structure (see also Martins et al. 12004(1 . We 
have tested the influence of those two effects in one of 
our low M models for HD 46202. Their combined effects 
lead to an increased ionisation in the outer atmosphere, 
the mean ionisation fraction of C IV being lowered by 
~ 0.1 dex (which does not modify the conclusions of 
Sect. I6.3|l . This slightly changes the UV line profiles, 
especially C iv AA1548,1551 which for a given M shows 
a smaller absorption in the bluest part of the profile. 
Quantitatively, the inclusion of advection and adiabatic 
cooling is equivalent to an increase of M by ~ 0.15 dex. 
We have thus included these two processes in our models 
for low density winds (HD 38666, HD 34078, HD 46202, 
HD 93028). 

Given the above discussions, we think our M deter- 
minations have a very conservative error bar of ±0.7 dex 
(or a factor ± 5). This is a quite large uncertainty which 
however does not modify qualitatively our results, namely 
the weakness of O dwarfs with low luminosity (see Sect. 



7. Discussion 

7.1. Evolutionary status 

Fig. 133 shows the HR diagram of the our sample stars. 
Overplotted on Fig. 1351 is our new calibration Tcs - lumi- 
nosity (Martins et al. 120051 solid line) for dwarfs: most 
stars of our sample agree more or less with this relation 
(within the error bars). The latest type stars of our sam- 
ple, which are also the stars showing the weakest winds, 
may be slightly younger than "standard" dwarfs of the 
same spectral type (or Toff). Notice that this docs not 
mean that these stars are the youngest in terms of ab- 
solute age, but that they are less evolved than classical 
dwarfs. Indeed, the youngest stars of our sample are those 



of Trumpler 16 (HD 93250 and HD 93204) for which we 
derive an age of 1 to 2 Myrs, compatible with the L— T^s 
relation. In comparison, HD 38666 and HD 34078 may be 

2 to 4 Myrs old according to our HR diagram (although 
given the error bars, we can not exclude younger ages), 
slightly less than for standard late type dwarfs. In the 
scenario where these two stars originated from a binary 
and were ejected in a dynamical interaction, Hoogerwerf 
et al. I|2()01(l estimate a travel time of 2.5 Myrs, while van 
Rensbergen et al. 1)1996(1 found travel times of ~ 3.5 Myrs 
for HD 38666 and ~ 2.5 Myrs for HD 34078. These es- 
timates are in good agreement with our results. We also 
derive an age of ~ 3-5 Myrs HD 46202, one of our weak 
wind stars. Note that this star is in the same cluster as 
HD 46223 which is likely 1-2 Myrs old according to Fig. 

03 A similar age should be expected for these two stars in 
case of a burst of star formation, but an age spread of 1-2 
Myrs (common in star clusters) can explain the difference. 
The same is true for the stars of Cr 228: HD 93146, the 
brightest star, may be slightly younger than HD 93028. 

HD 152590 behaves differently, being less luminous 
than other dwarfs of same Teff. It is interesting to note 
that this star is classified as Vz. Taking literally the re- 
sult of Fig. it seems indeed that it is younger than 
other dwarfs (but again, the error bars are large), con- 
firming the fact that Vz stars are supposed to lie closer 
to the ZAMS than typical dwarfs. However, HD 42088 is 
another Vz star of our sample, and it has a more standard 
position on the HR diagram. This poses the question of 
the exact evolutionary status of Vz stars. Indeed, they 
are defined by stars having He ii A4686 stronger than any 
other He II lines which is thought to be a characteristic of 
youth since this line is filled with wind emission when the 
star evolves. In fact the Vz characteristics may be more re- 
lated to the wind properties than to the youth of the star. 
Indeed, HD 42088 seems to have the same stellar prop- 
erties as HD 93146, but the former is classified Vz (not 
the latter) and has a weaker wind {M — 10^^ Mq yr^^ 
compared to lO"'^-^^ M© yr"! for HD 93146. Note how- 
ever that the distance (and thus luminosity) of HD 42088 
is highly uncertain. Obviously, more studies are required 
to better understand the physics of Vz stars. 

To summarise, there may be a hint of a link between a 
relative youth and the weakness of the wind if by youth we 
mean an evolutionary state earlier than for standard stars 
and not an absolute age, standard stars meaning stars with 
the average properties of dwarfs studied so far. But the 
present results are far from being conclusive. A forthcom- 
ing study of Vz stars in the LMC will probably shed more 
light on this issue. 

7.2. Wind properties 
7.2.1. Terminal velocities 

Some of the terminal velocities we derive are surprisingly 
low (see Table Ejl reaching values lower than the escape 
velocity in a one case (HD 34078). What could explain 
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Fig. 35. HR diagram of the Galactic stars. Evolutionary 
tracks are from Lejeune & Schaerer 1)2001(1 and Z = Zq. 
Isochrones for 0, 1.. 5 Myrs are indicated together with 
evolutionary paths of stars of different masses. The blue 
dashed line is our new relation Tcfi - luminosity for dwarfs 
(Martins et al. .2005,) . Squares are the stars studied here 
(open symbols are for Vz stars). 

this behaviour? First, the most obvious reason could be 
an underestimation of v^o ■ We have argued in paper I that 
in stars with weak winds the density in the upper parts 
of the atmosphere may be so low that almost no absorp- 
tion takes place in strong lines usually formed up to the 
top of the atmosphere. This explanation was also given by 
Howarth & Prinja H1989|l to justify the low Vq^ they ob- 
tained in some stars. If this is indeed the case, one would 
expect a smooth decrease of the absorption strength in 
the blue part of P-Cygni profiles due to the reduction of 
the density as we move outwards, and not a steep break 
as seen in dense winds. Is there such a transition? Fig. 
15^ shows the C iv AA1548,1551 line profiles of HD 34078 
and HD 46223 and reveals that although the increase of 
the flux level from the deepest absorption to the contin- 
uum level in the bluest part of the profile extends over a 
slightly larger range in the case of the weak wind star (3 A 
for HD 38666 instead of 2 A for HD 46223), it is difficult 
to draw any final conclusion as regards the reduction of 
the C IV AA1548,1551 absorption in the outer wind of low 
density wind stars from this simple eye estimation given 
also that blending is clearly apparent in the line of HD 
38666. More information is given by Fig. 1371 which shows 
the derived terminal velocities as a function of mean den- 
sity in the wind (see Sect. IH31 for definition). There is an 
obvious trend of lower terminal velocities with lower wind 
densities. This is not a proof of the fact that absorption 



Fig. 36. Comparison between C iv AA1548,1551 line pro- 
ffles in a star with weak (HD 38666, solid line) and strong 
(HD 46223, dotted line) wind. The rise of the flux level 
in the very blue part of the absorption features is slightly 
shallower in HD 38666. See text for discussion. 



in strong UV lines extends to larger velocities since low 
densities also mean low mass loss rate and correspond to 
stars with lower radiative acceleration. However, it is an 
indication that underestimations of Wqo are certainly more 
likely to happen in such low density stars. 

In view of the above discussion, it is not clear whether 
the lower density in the outer atmosphere of weak wind 
stars is responsible for an underestimation of the terminal 
velocities. But we can not exclude that our estimates of 
Vao are lower than the true values in low density winds. 
Now, with this in mind, let us now assume for a moment 
that the derived values are real terminal velocities: what 
are the implications? The radiation driven wind theory 
predicts that Voo is tightly correlated to the escape velocity 
(vesc) according to 

<^<cI{a)T^ (8) 

i — a 

where a is the usual parameter of the Castor, Abbott & 
Klein H1975|l formalism and I{a) is a correction factor to 
take into account effects of the finite cone angle of the star 
disk (see Kudritzki et al. I1989|l . In practice, it is possible 
to derive values of a from this equation once the stellar 
parameters and Wqo are known. The only problem comes 
from I{a) which is a complex function of (among other 
parameters) a. However, it is possible to solve this prob- 
lem with the following procedure: we first assume a given 
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Fig. 37. Terminal velocity as a function of mean density 
in the wind for the stars studied here. There is a clear 
trend of lower v^o in lower density winds. 

value of a, then estimate /(a) which is subsequently used 
to find a new a using 
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Fig. 38. Upper panel: ratio of terminal to escape veloc- 
ity in our sample stars (open symbols) and stars studied 
by Herrero et al. |200U1 1^001 Repolust et al. IMHI and 
Markova et al. 120041 ffilled symbols). The solid line indi- 
cates the classical value 2.6 derived for stars hotter than 
21000 K (Lamers et al. il995 ) . Dwarfs (giants, supergiants) 
are shown by triangles (squares, circles). Lower panel: de- 
rived a parameter from the estimated terminal and escape 
velocities. See text for discussion. 



A few iterations should lead to the final value of the a 
parameter. We have used such a scheme to estimate a for 
our sample stars and for a number of stars studied else- 
where (Herrero et al. .2000. .2002.. Repolust et al. 120041 and 
Markova et al. l2004|l . Solutions are usually found with less 
than 10 iterations (and are essentially independent on the 
starting value of a), except in the cases where Woo/vescwas 
larger than ~ 3; in that case, the iterative process did not 
converge but kept oscillating between two distinct values. 
The results for the cases where solutions could be found 
are shown in Fig. |2H1 (lower panel). A majority of cases 
lead to a ^ 0.5 — 0.6, in reasonable agreement with (al- 
though slightly lower than) theoretical expectations (Puis, 
Springmann & Lennon |2000|I . However, for the stars of 
this work with weak winds (and low Wqo), lower values 
are deduced {a 0.3). Hence, if the derived low terminal 
velocities correspond to real Voo, they may be due to low 
values of the a parameter (given Eq. |S1) If true, this may 
also have important implications for the scaling relations 
involving mass loss rates (see next section). Again, it may 
be possible that we underestimate the terminal velocities, 
but the above possibility is worth being discussed in view 
of the puzzle of the weak winds. Note that we have also 
plotted in Fig.|3Hlthe ratio of terminal to escape velocity 
which is usually of the order 2.6 for O stars with Tes > 



21000 K (Lamers et al. I1995f) . We see that hottest stars 
of our sample follow this general trend whereas stars with 
weak winds have much lower ratios. 

Another very interesting explanation for the low ter- 
minal velocities we derive (again assuming that they are 
the real values) is the effect of X-rays. Indeed, Drew et al. 
(1994) highlighted the fact the cooling time in the outer 
atmosphere of massive stars with relatively weak winds 
(late O / B stars) can become relatively high so that in 
the case where X-rays possibly emitted by shocks heat the 
outer atmosphere, this region remains hot. In that case, 
the ionisation structure is strongly modified compared to 
the inner atmosphere and in practice, the radiative force 
becomes negligible in this hot region. This means that the 
wind keeps expanding at the velocity reached at the top of 
the "cool" region which is lower than the value predicted 
by the radiation driven wind theory. This effect should be 
checked in future hydrodynamical simulations. 

7.2.2. Mass loss rates and modified wind momenta 

The mass loss of O stars has been known for a long time to 
depend on luminosity since due to the basic mechanism of 
radiatively driven winds, the more photons are available, 
the larger the acceleration and the larger the mass loss 
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rate (e.g. Castor, Abbott & Klein [WTSI Kudritzki & Puis 
I2000|l . Fig. 1211 shows mass loss rates for our sample stars 
(filled symbols) and stars from other studies (Herrero et al. 
l^nnnilM?^ Repolust et al. i2004i and Markova et al. 12001 
open symbols) as a function of luminosity. One sees that 
there is a good correlation between M and L for bright 
stars. Note however that our sample stars seem to show 
lower M than what could be expected from the other stud- 
ies (see also TableOl). For low luminosity stars, the correla- 
tion still exists, but the scatter is much larger. Moreover, 
the slope of the relation seems to be steeper for these 
objects, the transition luminosity being log ~ 5.2. 
Although our work is the first to show such a behaviour 
based on quantitative modelling of atmosphere of O stars, 
this trend was previously mentioned by Chlebowski & 
Garmany H19H1|I and Lamers & Cassinelli H199f)|l . This 
result confirms our finding of paper I in which we showed 
that the stellar components of the star forming region 
SMC-N81 displayed winds weaker than expected from the 
relation M - L at high luminosities. In paper I, we men- 
tioned that a possible explanation of such a weakness was 
the reduced metallicity of the SMC, but we also showed 
that the Galactic star 10 Lac had the same low mass loss 
rate. Here, we confirm that several Galactic stars with 
low luminosity indeed show low M, rendering unlikely the 
effect of metallicity alone. 

We also showed in paper I that the winds were weaker 
than predicted by the current hydrodynamical simula- 
tions. Fig. 0U| extends this trend for the Galactic stars 
studied here: in the "worst" cases, the difference between 
our derived M and the mass loss rates predicted by Vink 
et al. (j2()()l|l can reach 2 orders of magnitude! Note that 
even for bright stars our values are lower than the pre- 
dictions but only by a factor < 5. This is mainly due to 
the introduction of clumping in our models for these stars 
which naturally leads to reduced mass loss rates (Hillier 
et al.EOnS!- Bouret, Lanz & Hillier "SDITF). Using a different 
approach, Massa et al. (j2003,) also hilighted the possibility 
of lower mass loss rates due to the presence of clumping. 
We will come back to this below. 

One may also wonder why our values of M are lower 
than other previous studies. Indeed, as shown by Table 13 
our mass loss rates are systematically lower than derived 
so far for all stars of the sample. How can we explain 
this behaviour? First, let us recall that the mass loss rates 
gathered in Table |31 are estimated from either pure Ha 
analysis fLeitherer il98Hl Lamers & Leitherer ll998l Puis et 
al. [T996i. Repolust et al. l^Uni Markova et al. lMlHl or from 
pure UV analysis (Howarth & Prinia [r989l Chlebowski & 
Garmanv ll991|l . 

and Lamers & 



The Ha study of Leitherer (|T 
Leitherer H199;-{|l relied on measurement of Ha emission 
equivalent widths. They are linked with Ha luminosities 
which are themselves related to mass loss rates. The rela- 
tion L(Ha) - M is based on estimates of the population of 
the third level of Hydrogen for which departure coefficients 
from LTE are taken from the pure H He computations of 
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Fig. 39. Mass loss rates as a function of Luminosity for 
Galactic O stars. The filled triangles are the dwarfs stud- 
ied in the present paper (-1- 10 Lac from paper I displayed 
by the filled triangle without error bars). Open symbols 
are data from Herrero et al. ("2000" "2002), Repolust et al. 
I|2U(J4|I and Markova et al. (i2004,l . Triangles (squares, cir- 
cles) are for dwarfs (giants, supergiants) . 



Klein & Castor H1978|l . The Ha emission equivalent width 
is calculated from the total equivalent width to which a 
photospheric profile from the plane-parallel pure H He 
models of Auer & Mihalas H1972|l is subtracted. This pro- 
cedure may suffer from various approximations: the use of 
pure H He models may introduce errors in the prediction 
of departure coefficients since line-blanketing is known to 
affect the ionisation (and excitation) structure; moreover, 
photospheric profiles based on H He plane-parallel mod- 
els may also be different from line-blanketed spherically 
extended models. In the case of stars with Ha in ab- 
sorption, the estimate of the emission part may be risky 
since it may suffer from uncertainties in the photospheric 
component subtraction, from contamination by interstel- 
lar lines or from errors in normalisation process. Leitherer 
p988|l himself argues that the wind emission part of the 
global (wind -t- photospheric) profile becomes almost un- 
detectable in stars with Ha in absorption. Hence, the Ha 
mass loss rates of such objects based on this method is 
rather uncertain. 

Another method relying on Ha is that of Puis et al. 
I|1996(l and Markova et al. I|2()()4(l . It is again based on 
the emission part of Ha which is related to mass loss rate 
through an estimate of the H departure coefficient under 
the Sobolev approximation. Their method is accurate for 
values of their parameter A > 10^"* (see Eq. (3) of Puis 
et al. for a definition of A) which, for typical values of 
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Fig. 40. Comparison between derived mass loss rates and 
predictions of hydrodynamical simulations (Vink et al. 
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the stellar parameters of O stars, corresponds roughly to 
M > 10^^ M0 yr~^. This is mainly the reason why the 
authors give only upper limits for stars with Ha in ab- 
sorption. Finally, Repolust et al. (2004 ) used FASTWIND 
(see Sect. 16.2(1 to fit the Ha profile and estimate M. In 
the case of weak winds, this method becomes less and 
less accurate since Ha is almost entirely photospheric and 
hardly depend on mass loss. Here again, the authors only 
give upper limits on M in those cases. 

Concerning UV based determinations of the mass loss 
rate, Howarth & Prinja H1989|l use the column densities 
in several UV lines to estimate M x qi {qi being the ioni- 
sation fraction of the ion responsible for the line studied) . 
Under the approximation that the ionisation fraction is 
independent of luminosity, they derive M . The latter ap- 
proximation may introduce errors in the mass loss rates 
estimate. Indeed, our modelling of massive stars atmo- 
spheres reveals that the ionisation fractions are not con- 
stant when Tcft changes among dwarf stars (which in this 
case reduces to a change of luminosity, see Fig.|2J- Since 
the largest ionisation fractions of C IV given by Howarth 
& Prinja ((1989.1 are of the order 10~^-^ (see their Fig. 16) 
while we find values as high as ^ lO^"'^, a factor of 100 
between their M and ours is possible. 

The other UV analysis to which our results are com- 
pared is that of Chlebowski & Garmany H1991|l . The au- 
thors use fits of the UV lines using the method of Olson 
H1982|l which is similar to the SEI method (Lamers et 
al. I1987|l . Basically, this method uses a parameterisation 
of the optical depth through the line profile to produce 
synthetic profiles which, once compared with observed 



spectra, allow a determination of M. However, this lat- 
ter step requires a few approximations. In particular, the 
ionisation fraction has to be estimated which involves 
the use of SED at high energies (i.e. close to ionisation 
thresholds of C IV and N V): a blackbody distribution 
is used in the computations of Chlebowski & Garmany 
()1991(l . Moreover, only photoionisation and recombina- 



tions from/to the ground states are taken into account. 
Hence, once again the ionisation fractions may not be 
correctly predicted leading to errors on M. However, it 
is interesting to note that the approach of Chlebowski & 
Garmany ( l95T) is more accurate than that of Howarth & 
Prinja (^1989.) as regards the ionisation fractions and leads 
to lower mass loss rates (see Table , so that if, as we can 
expect, the ionisation fractions are better predicted in the 
current atmosphere models and are in fact higher, lower 
mass loss rates are not too surprising. 

Another important point highlighted in paper I was 
the behaviour of the so called modified wind momentum 
- luminosity relation (WLR) at low luminosities. Indeed, 
we showed that there seemed to be a breakdown of this 
relation below log ~ 5.2, at least for stars of the SMC 
(including the stars of paper I and 3 stars of NGC 346 
studied by Bouret et al. I2003|l . The Galactic star 10 Lac 
also showed a reduced wind momentum, indicating that 
this property could not be related to metallicity alone. In 
the present study, several characteristics of the WLR are 
highlighted. 

We first confirm that there is a breakdown - or at least 
a steepening- of the WLR below log ~ 5.2. Indeed, Fig. 
1421 shows that most stars below this transition luminosity 
have wind momenta lower that what one could expect 
from a simple extrapolation of the WLR for bright stars. 
Indeed, for log -r^ ^ 5.0, the relation for dwarfs -I- giants 



found by Repolust et al. (|2004|l gives wind momenta of 
the order lO'^^ while we find values as low as 10^^! There 
is in fact only one object which is marginally in agreement 
with the relation of Repolust et al. (I|2004() . but we have 
only an upper limit for M for this star (HD 152590). 

Second, we find that for the bright stars of this study, 
the modified wind momenta are reduced compared to the 
pure Ha analysis on which the WLR is established. The 
difference is on average a factor of between 5 and 10, 
especially for the two objects we have in common (HD 
15629 and HD 93250). The explanation of this discrep- 
ancy comes from the use of clumping in our models for 
these stars. Indeed, it is necessary to use inhomogeneous 
winds to correctly fit a number of UV lines, especially O V 
A1371 and N iv A1718. In Fig. EH we show the UV + Ha 
spectra of HD 93250 and two models: one with the mass 
loss rates derived by Repolust et al. I2004 II from Ha only 
(M = lO"^-'^'^ Mq yr~i), and the other with our esti- 
mate of M relying on both Ha and UV lines and taking 
clumping into account {M — 10~^'^^ Mq yr~^). One sees 
clearly that although both models are acceptable for Ha 
(in view of the nebular contamination one can not exclude 
one or the other possibility), UV lines are overpredicted 
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with M = 10-^ '^^ Mq yr^i. Note that Repolust et al. 
H2004|l put forward the fact that the presence of clumping 
may lead to an overestimation of the derived mass loss 
rates if unclumped models are used (see also Massa et al. 
I2(JU3|I . Their argument is mainly based on the larger modi- 
fied wind momenta of stars with Ha in emission compared 
to stars with Ha in absorption which can be explained by 
the neglect of clumping in the former. However, they do 
not exclude the existence of clumping in the latter, but 
claim that its effects can not be seen due to low optical 
depth. In our case, all stars have Ha in absorption, and 
we deduce the presence of clumping from UV lines. 

Note that in our study, we had to include clumping 
only in the more luminous stars to correctly fit the UV 
spectra. Does it mean that the winds of fainter stars are 
homogeneous? Not necessarily. Indeed, clumping is re- 
quired to reproduce O v A1371 and N iv A1718. But it 
turns out that in late type O stars (i.e. low L stars) O v 
A1371 is absent and N IV A1718 is mainly photospheric so 
that homogeneous winds give reasonable fits. Clumping 
may be present, but it can not be seen from the UV and 
optical spectra. In any case, if clumping were to be in- 
cluded in the models, the mass loss rates would have to be 
reduced to fit the observed spectra (see above). Thus, M 
would have to be further reduced compared to the already 
low values we obtain, making the winds of low luminosity 
O dwarfs even weaker! 

In spite of the global shift of our WLR for bright O 
dwarfs compared to pure Ha studies when clumping is in- 
cluded, the slope of the relation is roughly the same as that 
found by Repolust et al. I12()()4|I . This is important since it 
shows that the breakdown of the WLR we find at low lu- 
minosities is not an artifact of our method. Equivalently, 
this means that even if we underestimate the mass loss 
rates (due to ionisation fractions, see Sect. Vr>.[^ . there is 
however a qualitative change of the slope of the WLR near 
log 4^5.2. 

Can we estimate the value of the slope of the WLR 
in this low L range? The number of stars studied is still 
too low to give a reliable value, but if we do a simple by 
eye estimate, excluding star HD 152590 (due again to the 
fact that it is a possible member of a binary), we find a 
slope of the order of 4.3. As this slope is in fact equal to 
1/a' (e.g. Kudritzki & Puis .2000^1 . where a' = a — 5 and 
5 = 0.005. .0.1, we deduce a' ~ 0.25 and a - 0.30 which 
is very low compared to the classical value of ~ 0.6, but 
which is consistent with our finding based on the ratio of 
terminal to escape velocities fSect. 1772. l|l . This does not 
mean that a is indeed this low for these stars since both 
Voo and M may be underestimated, but it is at least a 
kind of consistency check. 

7.2.3. Origin of weak winds 

In view of the above results, what can we conclude as re- 
gards the origin of the weakness of the winds observed 
in some O dwarfs? The main possibilities have been de- 
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Fig. 41. Influence of clumping on the determination of M 
of HD 93250. The solid line is the observed spectrum, the 
dotted line is a model with M = 10~^'^^ yr~^(mass 
loss rate derived by Repolust et al. 120041 from Ha) and no 
clumping, and the dashed line is a model with M — 10~^'^^ 
Mq yr^^and /oo — 0.01. The inclusion of clumping leads 
to a reduction of the mass loss rate derived from both UV 
lines and Ha. See text for discussion. 



tailed in paper I. Among them, metallicity was the first to 
come to mind since at that time, as most stars with weak 
winds were found in the SMC (paper I, Bouret et dl.\^^. 
The present study clearly shows that metallicity cannot 
explain the reduced wind strength observed since several 
Galactic stars show mass loss rates and terminal veloci- 
ties as low as SMC objects. On the contrary, it becomes 
more and more evident that there is a transition in the 
wind properties near log = 5.2, although the reasons 
for such a change of behaviour remain unclear at present. 
To say things more clearly, we do not state that metal- 
licity has no effect on the wind strength (this is now well 
established): we simply show that low luminosity dwarf 
O stars have winds much weaker than expected from hy- 
drodynamical simulations and than so far observed for O 
stars, independent of metallicity. 

A possible explanation is the reduction of the a pa- 
rameter which, if it were as low as ~ 0.25, could explain 
both the reduced mass loss rates and terminal velocities. 
However, what could be the reason for such a low a ? 
A nice possibility was highlighted by Puis et al. H2000|l 
in their very detailed analysis of the line statistics. They 
first show that under fairly general conditions, the classi- 
cal a parameter, i.e. the one entering the slope of the line 
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Fig. 42. WLR for Galactic stars. Symbols and data are 
the same as in Fig. 13^1 Note the breakdown of the relation 
below log 4^ ^ 5.2. See text for discussion. The short 
(long) dashed line is the regression curve for supergiants 
(giants + dwarfs) of Repolust et al. H2004|l . 



strength distribution function, and the d parameter used 
to express the radiative acceleration according to 



Tad I - 

g oc t 



(10) 



in the CAK formahsm (see Castor, Abbott & Klein ll97^ 
are the same (this is the basics of the radiation driven wind 
formalism). However, the line strength distribution func- 
tion may not have a constant slope, and in that case the 
value of a such that a = a must be derived at the point 
where the line strength is equal to . This has two impor- 
tant consequences: first, since t = s^vth ^J/j^r (''^ith Se the 
Thomson scattering opacity and Vth the typical thermal 
velocity), it will be different from star to star; and second, 
t is not constant in a given star's atmosphere. In practise, 
this means that a should be different not only from star 
to star, but also from point to point in the atmosphere 
of a given star! However, Puis et al. (^20,00.) have shown 
that the slope of the line strength distribution function is 
constant over a large range of t values, implying that a 
keeps a constant value close to 0.6. But in extreme cases, 
we may reach a range where this slope varies: in that case, 
a is reduced. This is shown by Fig. 27 of Puis et al. ( 2000) 
where we see that for t lower than 10~^ a deviates signifi- 
cantly from 0.6. The interesting thing is that this situation 
corresponds to low densities (see the definition of t above) . 
In conclusion, a < 0.6 is expected in low density winds, 
which is consistent with our findings. 

Given this fact, the main problem would come from 
the strong disagreement between the results of spectro- 



scopic analysis and the predictions of hydrodynamical sim- 
ulations. Could such simulations overestimate mass loss 
rates? It is indeed possible since they neglect velocity cur- 
vature terms in the computation of radiative accelerations 
(Owocki & Puls ll99T?| . which can affect the final results. 
As discussed in paper I, O dwarfs with low luminosity are 
the most sensitive to such effects but test models for a 
40 M0 star performed by Owocki & Puis H1999|l lead to 
downward revision of the mass loss rate by only a factor 
1.5 while we would need a factor ~ 100! 

One can also wonder what is the effect of X-rays on 
the driving of winds. We have already seen in Sect. 17.2.11 
that low terminal velocities can be expected when X-rays 
are present due to changes in the outer atmosphere ioni- 
sation structure. Further insights can be found in Drew et 
al. H1994|l . These authors have studied two B giants and 
have found mass loss rates 5 times lower than values ex- 
pected for a simple extrapolation of the M - L relation 
of Garmany & Conti H1984|l established for O stars. They 
also detected X-ray emission in both stars, and argued 
that such X-rays, likely formed in the outer atmosphere, 
can propagate towards the inner atmosphere and change 
the ionisation structure here too, reducing the total ra- 
diative acceleration and thus the mass loss rate. This may 
partly explain why the mass loss predictions of Vink et 
al. I)2()()l|l . which does not take X-rays into account, are 
higher than our derived values. Obviously, hydrodynami- 
cal models including X-rays are needed to test this attrac- 
tive hypothesis. 

We highlighted in paper I that decoupling may be an 
alternative explanation although no conditions for it to 
take place were fulfilled in the N81 atmospheres. Here, 
the stellar and wind parameters of the weak winds stars 
being similar to those of the N81 stars, we have checked 
that such a process is not likely to be at work either. 

8. Conclusion 

We have derived the stellar and wind properties of 
Galactic O dwarfs with the aim of tracking the condi- 
tions under which weak winds such as observed in SMC- 
N81 (Martins et al. I2004|l develop. Atmosphere models 
including non-LTE treatment, spherical expansion and 
line-blanketing were computed with the code CMFGEN 
(Hillier & Miller ITMHIl . Optical H and He lines provided 
the stellar parameters while both UV lines and Ha were 
used to determine the wind properties. The main results 
can be summarised as follows: 

o The O dwarfs studied here are 1 to 2 Myrs old for the 
hottest and 2 to 4 Myrs old for the coolest. Except 
for the faintest, they have luminosities in reasonable 
agreement with the new calibration Tcff - Luminosity 
of Martins et al. ^UU^. 

o Stars with luminosities below a certain transition lu- 
minosity (log ^ 5.2) have mass loss rates of the or- 
der of l0~S- ~9-^ M0 yr~^ and low terminal velocities 
(down to 800 km s~^). The mass loss rates are lower by 
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nearly a factor of 100 compared to the hydrodynami- 
cal predictions of Vink et al. (2001). Uncertainties in 
the determination of M , discussed here in detail, are 
not expected to qualitatively alter the results. 

o Stars with log 4^ ^5.2 are found to have reduced 
mass loss rates compared to both hydrodynamical pre- 
dictions and previous analysis based only on Ha. The 
main reason is the inclusion of clumping in our models 
in order to fit O v A1371 and N iv A1718 in the lUE 
range. The adoption of pure Ha based mass loss rates 
does not allow fits of most of the UV lines. 

o The modified wind momentum - luminosity relation 
shows a break down around log -f^ = 5.2. Below this 
transition value, the slope corresponds to a value of 
the a parameter of the order of 0.3, which is consistent 
with the low terminal velocities observed. Such a low 
a is expected in low density winds (Puis et al. i2000J . 

o The origin of the weakness of the winds in low luminos- 
ity stars compared to hydrodynamical simulations is 
still unknown, but metallicity effects can be excluded 
since all the stars of the present study are Galactic 
stars and show reduced winds similar to SMC stars 
(Bouret et al. OIMI Martins et al. BDUljl . An earher 
evolutionary state than in standard dwarfs may be re- 
sponsible or not for the weakness: the present results 
can not resolve this issue given the error bars in the 
age estimates. 

o Although their origin remains unclear. X-rays appear 
to play a very important role in the physics of weak 
winds. They may be due to magnetic mechanisms and 
affect the ionisation structure. This can possibly re- 
duce the driving force and partly explain the low termi- 
nal velocities and low mass loss rates. Hydrodynamical 
simulations including X-rays should give more quanti- 
tative results. 

The low luminosity objects of our sample have not 
been studied individually with quantitative spectroscopy 
before since the atmosphere models allowing the analysis 
of weak wind stars have only been available for a few years. 
Indeed, the detailed modelling of UV wind sensitive lines 
requires a reliable treatment of line-blanketing since most 
of these lines are from elements heavier than He. This also 
explains why a number of previous quantitative analysis 
relied essentially on Ha. 

Now that the existence of weak winds has been estab- 
lished observationally both in the SMC and in the Galaxy, 
it would certainly be suitable to investigate the problem 
from a theoretical point of view with new hydrodynam- 
ical simulations. Apart from that, we still have to make 
sure that the ionisation fractions predicted by CMFGEN 
are correct since thay may alter the mass loss rate de- 
terminations. We will conduct such an investigation in a 
forthcoming paper using FUSE data. 
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